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A REVIEW OF SOME REYNOLDS NUMBER EFFECTS 
RELATED TO BODIES AT HIGH ANGLES OF ATTACK 

by 

Edward C. Polhamus 

Virginia Associated Research Campus? Newport News? VA. 


SUMMARY 

This paper presents a review of some effects of Reynolds number on 
selected aerodynamic characteristics of two- and three-dimensional bodies of 
various cross sections relavant to aircraft and missile fuselages at high 
angles of attack at subsonic and transonic speeds. Emphasis is placed on the 
Reynolds number ranges above the subcritical and angles of attack where lee- 
side vortex flow or unsteady wake type flows predominate. Lists of 
references, arranged in subject categories, are presented with emphasis on 
those which include data over a reasonable Reynolds number range, but also 
including some selected papers providing background information and 
prediction methods. Selected Reynolds number data representative of various 
aerodynamic flows around bodies are presented and analyzed and some effects 
of these flows on aerodynamic phenomena related to static longitudinal and 
directional chacteri sties, spin damping, and out-of-plane side forces are 
discussed. Finally, some suggestions for high Reynolds number research in 
the new National Transonic Facility on bodies at high angles of attack are 
offered. 
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INTRODUCTION 


The aerodynamic design of of modern aerospace vehicles is often 
critically dependent on the accuracy with which full scale viscous flow 
characteristics can be predicted. One of the areas that is particularly 
sensitive to Reynolds number effects is the high angle of attack 
characteristics of both aircraft and missiles. The current strong interest 
in providing expanded maneuver performance capbilities while maintaining 
satisfactory handling qualities and avoiding spin entry has increased the 
importance of high angle of attack aerodynamics and the need for full scale 
Reynolds number simulation. 

Current operational wind tunnels used for vehicle development fall short 
of simulating full scale Reynolds numbers by at least an order-of-magnitude 
for the many compressible flow conditions involved, thereby requiring large 
theoretical or empirical extrapolations which are often of questionable 
accuracy. This is particularly true for the high angle of attack conditions 
involving complicated separated flows with strong vortex interactions. This 
was one of the reasons for the decision to construct the National Transonic 
Facility (NTF) at the Langley Research Center. This high Reynolds number 
wind tunnel, based on the cryogenic pressure tunnel concept developed at 
Langley, is currently nearing operational status. 

As an aid to utilizing the capabilities of the National Transonic 
Facility in the most effective manner, it appears desirable to document the 
most pertinent Reynolds number studies available in the literature and 
analyze selected data of particular interest in relation to the anticipated 
research and development programs to be carried out in the NTF. Such 
documents should be of value in developing meaningful Reynolds number 
oriented research programs for the NTF, and thus make best use of the 
cryogenic pressurized tunnel concept, which can isolate certain parameters 
and provide new and improved simulation. 

The purpose of this paper, therefore, is to document one aspect of the 
high angle of attack aerodynamics described above; that of Reynolds number 
dependent flows related to the aerodynamics of bodies such as missiles and 
the fuselages of aircraft and spacecraft. 

The paper concentrates on the review and analysis of those data that 
cover a reasonable Reynolds number range and points out areas where 
extensions of such data using the capabilities of high Reynolds number 
2 



tunnels such as the NTF will be of particular value. Extensive reference 
lists are provided and divided into various categories as an aid to those 
wishing to make further use of the available literature on Reynolds number 
effects on bodies. While papers describing various prediction methods are 
listed}, no evaluation of their relative accuracies is attempted. Some papers 
which came to the author's attention after the text was completed are listed 
as "some additional papers." 

The discussion and reference lists are? in general, limited to the 
subsonic and transonic speed regimes. 

1. FORCES ON CYLINDERS IN NORMAL FLOW 

Many of the Reynolds number dependent viscous flow phenomena encountered 
by missiles and aircraft fuselages at high angles of attack can be related, 
at least qualitatively, to the flow about two dimensional cylinders in normal 
flow having the same or similar cross-sectional shapes. This qualitative 
similarity of the flows also extends into the compressible flow region. 

Because of this relationship between two-dimensional and three- 
dimensional flows about fuselage-type bodies and the fewer number of 
parameters involved in the two-dimensional case, it has been selected as the 
initial section of the present analysis. 

1.1 Circular cylinders 

In this section, the available data regarding the major effects of 
Reynolds number on the aerodynamic flow characteristics for circular 
cylinders immersed in a flow normal to the cylinder axis will be analyzed. 
Some of the general papers which provide good fundamental descriptions of the 
flow details are listed as references 1.1 through 1.5. Papers containing the 
experimental data to be utilized in this review are listed as references 1.6 
through 1.44. For convenience, the reference list is divided by subjects^ 
and within each subject they are, in general, listed in chronological order 
for historical interest. 

Since the emphasis in this paper is on the higher Reynolds numbers above 
the critical range, no attempt is made to cite the many low Reynolds number 
studies. For those interested in this range, Morkovin (reference 1.3) 
provides a rather extensive list of papers. 
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1.1.1 Smooth surface cylinders 


The various Reynolds number ranges to be discussed are identified in 
figure 1.1 and related to the separation, vortex shedding characteristics, 
and relative drag levels of smooth circular cylinders in incompressible 
flow. 

The "subcritical" range is defined as that range below a Reynolds 
number, based on cylinder diameter, of R D =2xl0 5 . In this range, laminar 
separation occurs on the front at a separation angle, 0 S , of approximately 

80° resulting in a wide wake and a high drag coefficient. In this range, the 
vortex shedding in the wake, while not of the classical "Karman" type, can be 
characterized as having a dominant periodicity. This range is followed by 
the "critical" range in which the laminar separation is followed by a 
turbulent reattachment enclosing a laminar bubble. The turbulent flow then 
separates further aft, at about 6=130°, at a Reynolds number of about 4xl0 5 , 
resulting in a narrowing of the wake and a large reduction in drag. In this 
range, the wake is unstructured with no dominant periodicity, and the vortex 
shedding frequencies become of the "wide-band random" type. 

The "supercritical" Reynolds number range, as defined here, extends from 
about R=4xlQ 5 to about R=6xl0 6 . In this range, transition moves upstream, 
forward of the laminar separation point and eventually to the vicinity of the 
forward stagnation point with the separation point moving forward from about 
130° to 115°, and there is a corresponding increase in drag coefficient. In 
the upper part of the supercritical region, the vortex shedding shifts to a 
"narrow-band random" pattern and the drag starts to level off. 

In the range above R=6xl0 6 , the vortex wake flow becomes at least 
"qua si-periodic" and some evidence of a drag reduction has been observed. In 
this paper, the term "hypercritical" is adopted for this range. 

It will be noted that the term "hypercritical" departs from historical 
precedence. The term "transcritical" has been accepted through usage for 
this region, although some controversy has appeared in the literature 
associated with its somewhat inconsistent connotation in describing the 
Reynolds number beyond the "supercritical" range. With the expected increase 
in high Reynolds number research in both two-dimensional and three- 
dimensional flows as a result of the availability of new high Reynolds number 
wind tunnels, consideration of adoption of a new term appears timely. In 
attempting to accomplish this in the past, some researchers have combined the 
4 



"critical " and "supercritical" regions and substituted the term 
“transcritical " followed then by "supercritical," While this provides a 
somewhat more logical terminology, the loss of the term "critical" describing 
the region where the separation characteristics are changing very rapidly and 
the shifting of the term "supercritical" from the region which it has defined 
for many decades appears undesirable. Therefore, since "transcritical" is a 
relatively new term whose use has been somewhat restricted due to the rather 
limited research in this region, this paper drops its use, substituting 
"hypercritical" and retains the well-established terms in the lower regimes. 
With the anticipated increase in two-dimensional and three-dimensional 
research in the "hypercritical" range, a consistent adoption of this term 
should soon eliminate the confusion that may result initially. 

Many correlations have been made of the effect of Reynolds number on the 
drag of circular cylinders in incompressible flow. However these 
correlations invariably have shown extremely wide bands of uncertainty, 
particularly in the higher Reynolds number ranges (see reference 1.5 for 
example). While the effects of surface condition, compressibility, end 
conditions, etc., are sometimes cited as possible causes; no concerted 
attempt has been made to isolate these effects, and therefore, such an 
attempt will be made here. 

The data selected from the literature search to establish the low speed 
drag coefficient for two-dimensional smooth circular cylinders is presented 
in figure 1.2 as a function of Reynolds number based on free stream 
conditions and the cylinder diameter. To ensure that effects of 
compressibility are insignificant, the data used has been limited to that 
obtained at Mach numbers less than 0,15. Comparisons of the drag variation 
with Reynolds number for cylinders of various diameters where the Reynolds 
number variations were obtained by velocity changes (references 1.9, 1.10, 
and 1.14) , indicated that measurable compressibility effects began in the 
range of Mach number from 0.18 to 0,20 with the value of 0.15 being selected 
as a conservative upper limit for esentially incompressible flow. 

To establish the level of drag coefficient at subcritical Reynolds 
numbers, data from the classical experiment performed by Wieselsberger at 
Gottingen (reference 1.6) is used. Wieselsberger tested 9 cylinders having 
diameters ranging from 0.05 to 300mm to cover the range of Reynolds numbers 
from 0.2 to 800,000. However, for the range of interest here only the data 
obtained with the 80mm (3.15 in) diameter cylinder is presented in figure 
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1.2. It will be noted that in addition to establishing the constant drag 
level for the subcritical range , his data establishes the beginning of the 
critical range* Although Wieselsberger also tested a 300mm cylinder to cover 
the full critical Reynolds number range, the data is not used here in view of 
possible end-plate boundary-layer and blockage effects associated with this 
extremely large diameter cylinder of relatively short length. 

To establish the remainder of the critical Reynolds number range, the 4" 
diameter cylinder data of Stack (reference 1.8), the 8.34” diameter of 
Pol ham us (reference 1.10) f and the 12" diameter of Del any and Sorensen 
(reference 1.9) are added to that of Wieselsberger. (It should be noted that 
in the subcritical range, the two data points from reference 1.10 are low 
due, in all probability, to inaccuracy associated with the very low dynamic 
pressures) . In the critical range, good agreement exists between all of the 
data, and it appears that the variation is accurately established. 

Although Stack obtained data in the upper critical range, the Mach 
number exceeded the 0.15 limit selected for the correlation. The Del any and 
Sorensen data in the lower critical range has been eliminated due to the 
effects of leakage at the wall pointed out by the authors. Their data in the 
upper critical and supercritical ranges is tentatively assumed valid in view 
of the greatly reduced end effects of three-dimensional cylinders in this 
Reynolds number range (see references 1.42 and 1.44). 

In the lower supercritical range, the data of Delaney and Sorensen, and 
Polhamus are presented in figure 1.2 and reasonably good agreement is 
observed. It should be noted that the data at higher supercritical Reynolds 
number presented in references 1.9 and 1.10 have been eliminated to avoid 
exceeding the Mach number limitation described earlier. The minimum drag 
coefficient, which occurs in the lew supercritical range, is about 0.25 and 
is in good agreement with the flight test data published by Stephens 
(reference 1.7). 

Data from other sources in this range have been omitted for the 
following reasons. The data of Humphreys (reference 1.11) is not utilized 
since the critical Reynolds number is low, indicating the probability of too 
high a stream turbulence level. In addition, only one data point below a 
Mach number of 0.15 and with the ends sealed was obtained in the 
supercritical range. The data of Jones et. al. (reference 1.15) in the lower 
supercritical range were obtained in freon at a tunnel stagnation pressure of 
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1/10 atmosphere and the resulting low dynamic pressure caused lew measurement 
accuracy (which was confirmed by the senior author). 

In the higher supercritical and hypercritical range of Reynolds number , 
it appears from the literature survey that only the data of Jones, et. al. 
(reference 1.15) satisfied the low Mach number and smooth surface condition 
requirements. Additional data sources often used in attempting to define the 
drag variation in the upper supercritical and hypercritical Reynolds number 
range are those of Roshko (reference 1.20) and Schmidt (reference 1.13). 
However, in the present study Roshko* s data, because of the surface condition 
of the cylinder, is included with the rough cylinder data rather than the 
smooth cylinder data. The Schmidt data was not utilized since it was 
obtained with a three-dimensional cylinder. Care must exercised in analyzing 
Schmidt's data due to the large variation in Mach number for a given Reynolds 
number. Results from the study by James et. al. (reference 1.27) did not 
meet the Mach number criteria. 

In the study by Jones et. al. , a 36" cylinder was tested in air for the 
upper supercritical Reynolds number range, and in freon to extend the data 
into the hypercritical range. These data, taken from table I of reference 
1.15, illustrates the increase in drag in the upper supercritical region as 
the transition point moves to the front face, resulting in a wider wake with 
less pressure recovery. The turbulent boundary layer can better withstand 
the adverse pressure gradient and the turbulent separation point remains aft 
of the laminar separation point. Above a Reynolds number of about 6x10 , in 
the so-called hypercritical Reynolds number range, the drag remains 
essentially constant as the transition point moves well forward on the front 
face of the cylinder and the separation point becomes relatively insensitive 
to Reynolds number. It is anticipated with further increases in Reynolds 
number that some decrease in drag will occur as fully turbulent flow is 
reached and the turbulent boundary layer begins to thin. A reduction of 
friction drag would be expected as well as an aft movement of the turbulent 
separation point. Tests in the new high Reynolds number wind tunnels would 
be useful in this regard. 

From the above analysis, it appears that much of the data scatter in the 
supercritical and hypercritical Reynolds number ranges shewn in previous 
summary papers can be identified as the effects of compressibility and 
surface condition, and that figure 1.2 represents a reasonably accurate drag 
variation for smooth circular cylinders in incompressible flew for Reynolds 
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numbers f rcm the subcritical through the supercritical range. As pointed 
out, additional data is required to define the variation in the hypercritical 
Reynolds number range, 

1.1.2 Effects of surface roughness 

The effect of surface roughness on the drag of circular cylinders normal 
to an incompressible flow is illustrated in figure 1.3. Since the primary 
interest in this paper is in the supercritical and hypercritical Reynolds 
number range, the data from the well-known early study carried out at lower 
Reynolds numbers by Fage and Warsap (reference 1.18) is not presented. In 
figure 1.3, the smooth cylinder "baseline" curve is taken from the 
correlation of figure 1.2, discussed in the previous section, in which the 
surface roughnesses selected appeared to be the lowest available and 
essentially "aer ©dynamically " smooth. It is recognized, however, that some 
small level of roughness drag may be present in the baseline curve, 
particularly in the hypercritical Reynolds number range. The roughness level 
for the various Reynolds number ranges is indicated for the baseline curve in 
terms of effective sand-grain height. 

The sources of data for rough cylinders at Mach numbers equal to or less 
than 0.15 which exceed the critical Reynolds number of the smooth cylinder 
are Roshko (reference 1.20), Achenbach (reference 1.21), Schulz and Hayn 
(reference 1.22), and Szechenyi (reference 1.24 and 1.25). 

The cylinder used by Roshko was an 18" diameter seamless black steel 
pipe which, after being cleaned by sandblasting, had a surface roughness of 
"about 200xl0“ 6 inches" resulting in a non-dimensional sand-grain roughness 
parameter, r/D, of about 10*“ 5 . Achenbach used glued spheres and applied a 
factor of 0.55 (obtained by curve fitting with the Fage and Warsap data in 
the critical range) to the sphere diameters to obtain equivalent sand-grain 
roughness (r s /D) values of 110xl0~ 5 , 450xl0 -5 , and 900xl0~ 5 . Szechenyi 1 s 
cylinders were roughened with calibrated spherical glass beads, and for the 
purpose of this paper, the factor of 0.55 used by Achenbach has been applied 
to convert Szechenyi* s data to an equivalent sand-grain roughness. It is 
recognized that such conversions are questionable and that the factor may not 
be independent of Reynolds number. However, since sand-grain roughness is 
more closely related to surfaces encounted on actual aircraft, the factor of 
0,55 is used as a "tentative" conversion, subject to further investigation, 
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In addition, the effects of tunnel turbulence variations, differences in 
roughness structure, etc,, makes the definition somewhat qualitative at 
best. However, it is thought that the data as presented illustrates the 
relationship between roughness and Reynolds number with regard to drag. 

Since the drag curves obtained by Achenbach are well defined by closely 
spaced data points, only his final curves are presented. In connection with 
the data from Szechenyi, it should be noted that some inaccuracy may have 
been encountered in transcribing and coverting the data from figure 3 of 
reference 1.25. 

Achenbach* s data is for relatively large roughness grain size and is in 
general good agreement with the early Fage and Warsap data with regard to the 
classical effect of roughness. This classical effect is to promote an 
earlier transition to a turbulent boundary layer, and thus a reduced pressure 
drag. This reduction is less than that obtained on smooth cylinders in the 
critical range and is followed by a more rapid increase in the supercritical 
range. It has been suggested by Batham (reference 1.23) that the higher 
value of minimum drag is due to the roughness dispersing the turbulent bursts 
and leaving the boundary layer less capable of withstanding the adverse 
pressure gradient. Achenbach* s study extends the data on highly rough 
cylinders into the upper supercritical Reynolds number range and illustrates 
the leveling off of the drag in this region as a value approaching that for 
the subcritical smooth cylinder. This is the result of the transition point 
moving close to the stagnation point so that separation occurs near the 
minimum pressure point similar to the subcritical laminar separation case. 

For large values of roughness, the drag becomes essentially indepedent of 

r/D. 

For smaller size roughness, the data of Szechenyi appears to indicate 

that below the supercritical Reynolds number range, roughness levels less 

than r/D=8xlO - ^ probably do not produce any measurable effect on the drag. 

As the boundary layer thins in the supercritical range, however, the 

roughness level of r/D=8xlO” 5 begins to produce a measurable drag increment 

for R d values greater than about 10 . This trend appears to be substantiated 

by Roshko*s data in the upper supercritical and hypercritical range where an 

appreciable drag increment due to roughness is observed for a roughness 

— s 

parameter of only 1.1x10 
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The dashed lines used in figure 1.3 to connect the low roughness data to 
the smooth cylinder curve are entirely conjectural but are presented to 
illustrate the possible onset of roughness effects above the critical 
Reynolds number. 

With regard to the drag of rough cylinders in the hypercritical Reynolds 
number range, it is of interest that data was obtained in this high Reynolds 
number range as early as 1929 by Dryden and Hill and published in a seldom 
cited paper (reference 1.19). They obtained data, in natural winds, on a 
power plant chimney constructed of brick. But due to three-dimensional 
effects and undetermined equivalent sand type roughness, no attempt is made 
here to include the data in figure 1.3. 

From the results of figure 1.3, it appears desirable to establish the 
onset of roughness effects in the Reynolds number range above 10^. This may 
require a smooth reference test with surface finish of higher polish than 
tested previously. 

Additional data on the effects of surface conditions in compressible 
flow would also be of value particularly in the subsonic range and at 
Reynolds numbers above 10®, Data obtained at lower Reynolds numbers suggest 
that for compressible flew at subcritical Mach numbers, surface roughness may 
reduce the drag even in the supercritical Reynolds number range (see 
references 1,22 and 1.27). For supercritical Mach numbers, with shock- 
induced separation on the forward face, Mach a (references 1.36 and 1,37) has 
indicated drag reductions due to surface roughness for Mach numbers up to 0,9 
over the Reynolds number range for which data is available. 

1.1.3 Effects of compressibility 

As mentioned earlier, measurable effects of compressibility are 
encountered at Mach numbers well below the critical Mach number which? 
depending upon the boundary layer condition, occurs in the range of 0.40. 

This is not surprising since at a Mach number as low as 0.25, the minimum 
pressure coefficient is more negative than the incompressible value by nearly 
9% according to the Karman-Tsien equation (reference 1.30). This increases 
the magnitude of the adverse pressure gradient thereby moving the separation 
point forward and reducing the pressure in the separated region. Figure 1.4 
summarizes the existing smooth cylinder data for Mach numbers in the range 
from essentially incompressible to 0.50 which is slightly in excess of the 
10 



critical. For applications related to the subject of this paper , Mach 
numbers normal to the cylinder axis up to 0.5 are usually sufficient. If 
information for higher Mach numbers is required, data can be found in 
references 1.33, 1.34, 1.36, and 1.37, although the data is limited to 
Reynolds numbers below 10 . 

In figure 1.4, the line representing the incompressible variation of 
circular cylinder drag coefficient with Reynolds number is from figure 1.2 
and represents Mach numbers equal to or less than 0.15. In the subcritical 
Reynolds number range where the incompressible laminar separation point 
occurs on the front face, the effects of compressibility for Mach numbers up 
to 0.5 are relatively small and appear to be well established by the data of 
Stack (reference 1.31) and Gowen and Perkins (reference 1,34). However, in 
the higher Reynolds number range of interest for this paper, the trends with 
Mach number can only be partially established. 

From the data obtained by Stack (reference 1.31) and Schulz and Hayn 
(reference 1.14) on cylinders of various diameters, it is fairly well 
established that at least below the critical Mach number, compressibility 
shifts the critical Reynolds number range to considerably higher Reynolds 
numbers as seen in figure 1.4. It is presumed that this effect is associated 
with a delay in the formation of boundary layer turbulence due to the 
increased magnitude of the favorable pressure gradient on the forward face, 
and an increase in the adverse gradient aft of the minimum pressure point. 
That the changes in pressure gradients can be appreciable is seen in the fact 
that the minimum pressure coefficient is more negative by nearly 24% at a 
Mach number of 0.40 according to calculations by the inviscid full-potential 
transonic theory described in reference 1.39. These changes in pressure 
gradient affect not only the separation point and wake width but also result 
in a more negative pressure at the separation point, and therefore, in the 
wake. This results in an increased drag over the entire Reynolds number 
range. 

The only applicable data on compressibility effects in the upper 
supercritical and lower hypercritical Reynolds number range appears to be 
that obtained by Jones et. al. (Table II of reference 1.35) , and their data 
is shown in figure 1.4 for constant Mach numbers of 0.25, 0,35, and 0.40. 

For the Mach number of 0,25 both the air and freon data is presented, but for 
the other two Mach numbers the freon data has teen omitted in view of 
possible real gas effects which tend to cause departures from air results in 
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the higher local Mach number conditions (see references 1,40 and 1.41), In 
this regard*, it is of interest to note that the rather unusual peaky type 
variation of c^ with Reynolds number for M=0 .35 and M=0„4Q shown in figure 12 
of reference 1.35 would essentially disappear if the air and freon data were 
separated and if the 1/10 atmosphere freon data (with its low accuracy) was 
eliminated. Since the freon data represents the highest Reynolds number data 
available, the degree to which it represents results in air might be worthy 
of further study. 

Although the data is rather sparce, it appears from figure 1.4 that the 
effect of compressibility at subcritical Mach numbers is a minimum at low 
supercritical Reynolds numbers and increases considerably in the 
hypercritical range as the boundary layer becomes thinner. 

Based on the data of figure 1.4 it is apparent, that except for the data 
of reference 1.35 in the lower hypercritical range, very little penetration 
has been made into the Mach-Reynolds number envelope above a Reynolds number 
of about 5x10^ , The new cryogenic wind tunnels provide a source for the 
needed data. 


1.1.4 Effects of fineness ratio 

Although fineness ratio ( length- to-diameter ratio, L/D) effects involve 
three-dimensional flows, the cylinders (at a=90°) with "squared-off" ends 
will be covered here while the conventional nose shapes appropriate to 
fuselages will be treated in a subsequent section on three-dimensional 
bodies. 

The flow around the ends of cylinders of finite length has a pronounced 
effect on the drag. Possible sources of the drag changes include various 
combinations of venting from the end region into the near wake (which 
increases the pressure in the separated region) reduced velocities and 
Reynolds numbers near the ends, and separation at the sharp corners of the 
"squared-off" tips. The primary sources of data for the effect of fineness 
ratio are the experimental investigations of Wieselsberger (reference 1.42) 
and McKinney (reference 1.44). Both investigations show a very large 
reduction in drag, at subcritical Reynolds numbers, as the fineness ratio is 
reduced which is attributed primarily to the aforementioned venting into the 
near wake. 
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To investigate the influence of Reynolds number? Wieselsberger tested 
one finite length cylinder through the critical range* and McKinney 
investigated a large number of fineness ratios over a wide range of Reynolds 
numbers. Reynolds number variations for two fineness ratios selected from 
these studies are shown in figure 1.5 in terms of the drag proportionality 
factor? n , (ratio of finite length cylinder drag coefficient to that for an 
infinite length cylinder) in the critical Reynolds number range (the flagged 
symbol is from Wieselsberger and is used to extrapolate McKinney's data to 
the subcritical level) . The results indicate that increases with Reynolds 
number and approaches a value of 1.0 at the end of the critical range. This 
is apparently associated with a reduction of venting into the near wake for 
the low drag conditions where the wake is narrow and the pressure recovery 
rather high. McKinney has shown that under these conditions, lower fineness 
ratios can generate drag coefficients which are higher than the two- 
dimensional case (n>1.0) due, at least in part, to the increased importance 
of the flow separation around the "squared-off" ends (see reference 1.44). 

That portion of McKinney's supercritical data which was obtained at Mach 
numbers less than 0.15 is also shown in figure 1.5 and illustrates the 
reduction of n as the turbulent separation point moves forward widening the 
wake and reducing the pressure recovery, and thereby again allowing 
pressurization due to end effects. 

The only data found in the upper supercritical Reynolds number range is 
that obtained by Dryden and Hill on a 10-foot diameter cylinder in natural 
winds (reference 1.19). The factor for this data was determined by comparing 
the three-dimensional data extracted from the natural wind data with the two- 
dimensional correlation from wind tunnel data presented in figure 1.2. 

Because of uncertainties regarding the turbulence level, velocity 
distribution and surface condition for the natural wind tests the extracted 
factor is shown as a band. Despite these uncertainties, the values of 
are such that they reinforce the expectation that drag reductions due to 
finite length continue into the upper supercritical range. 

The reader is referred to reference 1 .44 for the effects of Reynolds 
number on the drag of finite length cylinders with hemispherical ends or 
noncircular cross sections. 
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1.1,5 Lifting cylinders 


An investigation of the effects of Reynolds number on the aerodynamic 
characteristics of two-dimensional lifting cylinders of circular cross 
section in incompressible flow has been carried out by Lockwood and McKinney 
(reference 1.11) and some of their results are shown in figure 1.6. This 
lift was generated by a spoiler type flap, and both the flap projection, h, 
and orientation angle, 6 , were investigated. For the data shown, the flap 
projection was 6% of the cylinder diameter. 

In figure 1.6, the lift coefficient, c^ , is presented as a function of 
Reynolds number, R D , for various flap orientation angles, 6 » The flap acts 
as a spoiler and tends to fix the separation line on the lower surface of the 
cylinder. The resulting wake deflection and lift force depends on the 
location of the upper surface separaton line. The results shown in figure 
1.6 illustrate the large increase in lift coefficient in the critical 
Reynolds number range for 5=90° from the free stream velocity vector) 
associated with the rearward movement of the turbulent separation line on the 
side opposite the flap resulting in large negative pressures which are 
augmented by the induced circulation. It will be noted that the maximum lift 
coefficient, which occurs in the upper critical Reynolds number range, 
approaches a value of 1.6, which is about 2.5 times the drag coefficient. 

As the Reynolds number increases in the supercritical range, the lift 
decreases due to the forward movement of the turbulent separation line. The 
loss of lift resulting generally from positioning the flap either fore or aft 
of the 90° position is also shown. For detailed pressure measurements and 
other flap projections, the reader is referred to reference 1.11. 

Although the above study was made to investigate the application of the 
lifting cylinder concept for the recovery of expended cylindrical rocket 
boosters, the results are also related to the use of slender nose strakes on 
three-dimensional fuselages at moderate and high angles of attack. For 
example, similar Reynolds number dependency might also be expected with 
retractable nose strakes used to develop side forces at high angles of attack 
to augment declining directional control. The potential hazard of developing 
such controls in the wind tunnel at subcritical Reynolds numbers is 
apparent. 

A somewhat more familiar application of slender nose strakes is that of 
providing directional stability at high angles of attack to offset the loss 
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of vertical tail effectiveness. In this application, the high crossflow 
velocity at the windward strake relative to that on the leeward strake 
results in an increased vortex strength and a net vortex induced side force 
that augments the directional stability. The magnitude of the forces 
generated can be influenced by Reynolds number through its effect on the 
natural separation location and the velocity at the strake location. 

(Limited data on two-dimensional cylinders with two flaps 180° apart that is 
related to the above strake application is also contained in the Lockwood and 
McKi nney st udy . ) 

1.2 Non-Circular Cylinders 

For a variety of practical reasons, the cross-sectional shapes of 
aircraft fuselages often deviate from a circle. For modern combat aircraft 
with relatively long fuselage forebodies ahead of the wing and having the 
requirement to maneuver to high angles of attack, the crossflow 
characteristics of their non-circulator forebodies often play an important 
role in the aerodynamic characteristics of the aircraft. This is 
particularly true with regard to such handling qualities as directional 
stability, nose slice, spin damping, pitch-up, etc. 

Shortly after the development of Allen' s crossflow theory (see section 
4) , it became apparent that data on two-dimensional cylinders of various non- 
circular cross-sections was needed in order to apply the method to non- 
circular shapes. Although some data was available, it was generally limited 
to subcritical Reynolds numbers. This prompted the two-dimensional cylinder 
drag study carried out by Del any and Sorensen (reference 1.45) that covered a 
variety of shapes over a range of Reynolds numbers. Their tests were made at 
zero flow incidence angle, and therefore, limited primarily to applications 
related to the longitudinal aerodynamics of three-dimensional bodies. 

Later, a low speed study by Letko and Williams (reference 1.46) 
illustrated the strong effect of fuselage cross-sectional shape on the 
directional stability of aircraft type configurations at moderate and high 
angles of attack. The study of cross-sectional shape on aircraft stability 
was soon extended to additional shapes and to higher speeds by other 
researchers (see references 1,47 through 1.50). The large effects 
encountered and the fact that the separated flows involved would be Reynolds 
number dependent prompted the two-dimensional non-circular cylinder studies 
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by Polhamus et. al. (reference 1.51 and 1.52) and Lockwood (reference 1.53) 
in which the effect of flow incidence was varied and the forces parallel and 
perpendicular to the wind were measured. These data, which covered 
subcritical, critical, and supercritical Reynolds number conditions, allowed 
application of the crossflow theory for the tested geometries to combined 
angles of attack and sideslip and to combinations of spin rates and sink 
speeds. 

The purpose of this section on non-circular cylinders is to present some 
representative examples of the important Reynolds number effects encountered 
along with a brief discussion of the implication, particularly with regard to 
future aerodynamic research. The reader is referred to references 1.45 and 
1.51 through 1.53 for the complete data presentations and discussions. 

1.2.1 Zero flow incidence angle 

The variation of drag coefficient with Reynolds number for some selected 
two-dimensional cylinders of various non-circular cross-sectional shapes at 
zero flow incidence is presented in figure 1.7 through 1.10. Since all of 
the shapes presented are flat on the windward face, the Reynolds number has 
been based on the maximum width of the cylinder normal to the free stream 
direction since it is a somewhat more meaningful length for these shapes with 
regard to boundary layer development as the windward corners are approached. 

Figure 1.7 summarizes the low subsonic speed data from references 1.45, 
1.51, and 1.52 for square cross-sections having various values of the non- 
dimensional corner radius, r/w, where r is the radius and w is the width of 
the cylinder. The results illustrate the large effect of corner radius both 
on the level of the drag and on the Reynolds number for transition from 
subcritical to supercritical flow. For the smallest radius, r/w=.021, 
separation occurs on the windward corners, due to the strong adverse 
gradient, and the drag is essentially the same as that for a flat plate. 

Some drag reduction might be expected, even for this small radius, at some 
Reynolds number beyond the limit of the test data, and the extrapolation 
shown is based on the correlation presented in figure 9 of reference 1.51, 
which represents a conjectured trend. The reduction of drag in the 
subcritical Reynolds number range, as the radii are increased, appears to be 
associated primarily with some recovery of the suction pressure loops on the 
windward corners due to the reduced adverse pressure gradients. The earlier 
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occurrence of the critical Reynolds number range and the reduction in 
supercritical drag level also is associated with the reduction in adverse 
pressure gradients. A more detailed analysis can be made with the aid of the 
measured pressures for an r/w of 0.245 presented in Table I of reference 1.52 
and by comparisons with inviscid theoretical results, such as presented in 
reference 1.54, 

The fact that the drag variation for the square cylinder does not 
monotonically approach that for the circular cylinder as the corner radius 
increases is due to the fact that the square has four suction pressure loops 
compared to two for the circle. This and the fact that the drag of the 
diamond-shaped cylinders do approach the circle monotonically is discussed in 
reference 1.51. For a summary of the effects of Reynolds number on a series 
of diamond-shaped cylinders having various corner radii, the reader is 
referred to figure 9(b) of reference 1.51. 

An example of the effect of compressibility on the drag coefficient of a 
square cylinder is presented in figure 1.8, which was obtained from some 
unpublished measurements made in the Langley High Speed 7-by-10 Foot Tunnel 
on a semispan model having an effective fineness ratio of 16 and an r/w of 
0.245. 

Figure 1.9 illustrates the drag reduction, for given corner radii, 
associated with increases in the cylinder depth in the stream direction. 

This drag reduction is due primarily to improved flow reattachment downstream 
of the windward corners. 

The effect of corner radius on basic triangular cross-sections is 
presented in figure 1.10. In order to provide the probable lower bound over 
a wide Reynolds number range, the case for an r/w value of .500 (circle) is 
reproduced from figure 1.2. 

Although a considerable amount of data exists for non-circular cylinders 
at zero flow incidence, it is primarily limited to incompressible flow and 
Reynolds numbers less than 2x10 . Research needs to be extended into the 
compressible flow regime and to higher Reynolds numbers. 

1.2.2 Effect of flow incidence 

In references 1.51 through 1.53, non-circular cylinders were tested over 
a wide range of flow incidence angles and Reynolds numbers to identify the 
types of viscous dependent crossflows that might be expected on non-circular 
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forebodies at high angles of attack in sideslip* In general? the results 
indicated that large changes in cross wind or side force coefficient? Cy* 
are encountered that can have an important bearing on aircraft directional 
stability at high angles of attack and on spin damping. 

Some examples of the effects of cross-sectional shape obtained are 
presented in figures 1.11 through 1.14, where the side force coefficient is 
shown as a function of Reynolds number, R^ for a constant flow incidence 

angle, <j> . A value of <j> of 10° has been selected for these examples, and the 
reader is referred to references 1.51 through 1.53 for the results of other 
flow incidence angles and other cross-sectional shapes. 

The critical effect that Reynolds number can have on the side force 
generated by square cylinders with round corners is illustrated in figure 
1.11. Side force data for corner radii of 8% and 24.5% of the cylinder width 
are presented along with pressure distribution data for the 24.5% radius 
cylinder. Of course, for pure inviscid flow with no "Kutta condition" 
imposed, no side force would be developed. For the viscous flow case, 
however, large side forces can be developed which generally change sign with 
Reynolds number, as illustrated by the data for the 24.5% radius cylinder. In 
the lower Reynolds number range, a postive (towards the wind) side force is 
developed followed by a transition to a negative side force beginning at a 
Reynolds number, R^, of about 0.5x10®. The flow mechanism generating these 
forces can be seen from the subcritical and supercritical pressure 
distributions shown in the sketches. These pressure distributions were taken 
from the experimental results of reference 1.52 and are compared with the 
inviscid theory of reference 1.54. For the case shown, the terms 
"subcritical" and "supercritical" refer only to the flow conditions around 
the front leeward corner (left side in the sketch) . The pressure 
distribution obtained at a Reynolds number of 0.3x10® indicates laminar 
separation at the front leeward corner due to the strong adverse pressure 
gradient, with only a fraction of the theoretical suction loop being 
achieved? attached flow being developed around a sizable portion of the 
front windward corner, as a result of the reduced adverse gradient? and 
turbulent separation on the rear windward corner. The positive side force 
results from the combination of the near-theoretical suction pressures on the 
front windward corner, the partial suction pressure development around the 
rear windward corner, and the loss in suction pressures on the leeward 
corners. As the Reynolds number increases, the pressure loop development 
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essentially reverses as attached flow develops on the front leeward corner , 
with turbulent separation on the rear leeward corner. In addition* a 
circulation is established which reduces the suction pressure loop on the 
front windward corner and increases the magnitude of the front leeward 
suction pressure loop beyond that for the inviscid theory without 
circulation. As a result of these flow changes? a large negative side force 
is developed. This force is nearly equal to 2fTsin<l>? the normal force 
generated by a flat plate in potential flow at a flow incidence angle? ? and 
a Kutta condition at the trailing edge. 

The results shewn for the smaller (8%) radius square cylinder shown in 
figure 1.11 illustrate the higher level of Reynolds number required to 
produce attached flow around the front leeward corner as a result of the 
steeper adverse gradient associated with the smaller corner radius. 

The significance of the forces developed with regard to fuselage 
forebodies are that a positive side force translates into a restoring moment 
with regard to directional stability at high angles of attack and translates 
into a propelling moment during a flat spin. With the side force changing 
sign at supercritical Reynolds numbers? it is very important that full-scale 
Reynolds number conditions be accurately simulated when investigating high 
angle of attack directional stability and spin characteristics. This was 
discussed in detail in reference 1.51 and led to the current procedures used 
in spin tunnel testing. A more complete discussion of the effects of 
Reynolds number on directional stability and spin damping will be presented 
in section 4.2. 

According to figure 1.11 (using the 24.5% radius cylinder as an 
example) ? additional critical Reynolds numbers might be expected both below 
and above that shown. For example? at lower Reynolds numbers separation 
would be expected to occur well forward on both front corners and a side 
force value closer to zero would probably occur. At Reynolds numbers above 
the highest shown in figure 1.11? the suction loops around the rear corners 
might be expected to develop? first on the windward corner and then on the 
leeward corner. With these flow changes and the resulting circulation 
changes? additional variations in side force might be expected to occur. Such 
multiple critical Reynolds number effects would have important implications 
in relation to directional stability and spin damping. 
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The capabilities ot the new high Reynolds number cryogenic wind tunnels 
offer a means of exploring these passible multiple critical Reynolds number 
effects? and this will be discussed more in section 4, 

While a complete summary of the available two-dimensional data on non- 
circular cross sections is beyond the scope of this paper? a few more 
examples of the effects of Reynolds number on side force are presented in 
figures 1.12 through 1.14. Figure 1.12 illustrates the difference in force 
characteristics generated by a rectangular cross section having rounded 
corners and a cross section with a similar shape on the windward half but an 
elliptical shape on the leeward half (referred to as a triangular shape). 

The rectangular cylinder characteristics are similar to those of the 
previously discussed square cylinder. However? because of its steeper 
adverse pressure gradients? the triangular cylinder encounters the critical 
Reynolds number on both corners within the test range and at higher Reynolds 
numbers than for the rectangle. No additional large Reynolds number effects 
would be expected for the triangular cylinder, while additional effects 
smilar to those discussed relative to the square cylinder would be expected 
for the rectangular cylinder. 

Figure 1.13 illustrates the effect of corner radius for two basically 
triangular type cross sections taken from the results of references 1.51 and 
1.53. The shift of the positive side force loop to a higher Reynolds number 
for the cylinder having the smaller corner radii is due? of course? to the 
more severe adverse pressure gradients. 

An illustration of the critical nature of slight departures from a 
circular cross section on the side force developed in the critical number 
range is presented in figure 1.14. This data was taken from the study 
carried out by Lockwood (reference 1.53) on a parametric series which 
transitioned between the circular cylinder and cylinder E of figure 1.12. 

For the complete series results over a wide range of flow incidence angles? 
the reader is referred to reference 1.53. 
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FORCES ON CIRCULAR CYLINDERS IN OBLIQUE FLCM 


2 « 


Inasmuch as the flow at angle of attack or sideslip about a three- 
dimensional body of revolution having a cylindrical afterbody of circular 
cross section is related in several respects to that about a cylinder of 
infinite length in oblique flow? a brief review of the infinite cylinder 
characteristics is presented in the following sections. 

The reader will note that in the interest of historical clarity? the 
description of the early studies of oblique cylinders will maintain the 
references to sweep and yaw of interest at that time. However? in the 
discussions of applications more directly related to the subject of this 
paper? the angle of attack? a, is used to define the oblique angle. 

2.1 Independence Principle 

Spurred by the growing interest in swept wings as a means of delaying 
the onset of compressibility effects? theoretical studies of the boundary 
layer development on yawed plates and infinite cylinders were begun during 
the mid 1940 8 s as an initial step towards the development of a completely 
three-dimensional viscous flow theory. Studies of laminar boundary layers by 
Prandtl (reference 2.2) and Struminsky (reference 2.8)? using the Navier- 
Stokes equations? illustrated that the boundary layer characteristics in 
attached laminar flow are independent of the spanwise velocity component. 

This finding? coupled with the earlier discovery that the inviscid flow 
velocity perturbations due to thickness are also unaffected by the spanwise 
velocity component? set the stage for the important research that followed. 

Jones (reference 2.3)? in discussing what Prandtl 8 s work meant in terms 
of the boundary layer separation on yawed cylinders? coined the term 
"independence principle." Jones concluded that? for the laminar flow case at 
least? the flow around the cylinder should depend on conditions in the cross- 
flow plane only. In particular? he concluded that (a) the distance to the 
separation point in the crossflow plane should be independent of yaw angle? 

(b) the critical Reynolds number based on cylinder diameter and crossflow 
velocity should also not change with yaw angle (however? the critical 
Reynolds number based on cylinder diameter and freest ream velocity should 
increase? by 1/cos 8? with increasing yaw angle) ; (c) the "drag" force normal 
to the cylinder axis should be reduced by cos^8 as the yaw angle increased? 
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since the dynamic pressure in the crossflow plane would be reduced by that 
amount. Essentially,, these same conclusions were reached independently by 
Sears (reference 2.4) , Other early researchers in this area were Wild 
(reference 2.5) and Cooke (reference 2.6). 

2.2 Experimental Results 

The later experimental studies with yawed cylinders indicated that, 
while the independence principle generally applied in the low Reynolds number 
range where there was laminar separation which is relatively independent of 
Reynolds number, it did not apply at the higher Reynolds numbers where there 
was turbulent separation. The most complete experimental study of flow 
around yawed circular cylinders was that by Bursnall and Loftin (reference 
2.7) , which covered both critical and supercritical flow conditions. Their 
study showed that the critical Reynolds number, rather than being delayed to 
a higher freestream velocity as the independence principle indicated it 
should, actually occurred at lower freestream velocities as the yaw angle 
increased. From this study and that of smith et. al. (reference 2.13), it 
appears that the transition from laminar to turbulent flow in the near wake 
and forward to the laminar separation point is indeed affected by the axial 
flow component. This results from the effect of the axial component in 
diverting the streamlines to a path more parallel to the freestream direction 
than normal to the axis. With the boundary state being related more to the 
"effective" Reynolds number along the streamline, Bursnall and Loftin 
indicated that the critical Reyolds number is nearly independent of yaw angle 
if it is based on chord lengths and velocities in the freestream direction 
rather than on those in the crossflow direction. Several methods have been 
proposed for approximating the "effective" Reynolds number related to the 
streamline path (see references 2.14 through 2.16, for example)? among these 
it appears that Esch® s factor (reference 2.15) collapses the Bursnall and 
Loftin data somewhat better than do the others and will be reviewed below. 

\ 

In the following, the oblique angle will be defined as angle of attack, a, in 
keeping with the subject of this paper. 

The Esch factor approximates the angle- of- attack effects on streamline 
length by use of the ellipse formed by the intersection of the surface and 
the plane parallel to the freestream that cuts through the cylinder at an 
angle, a , with respect to its axis. By use of this concept, an "effective" 
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Reynolds number, Rgjf, can be related to that based on diameter, R D , as 
follows; 


R eff R D R ot 


where K a = 


0.983/sina 


+ 0.311 + 0.287 sina 


1.581 


In the above formulation, Ka is equal to the ratio of the 
circumferential length in the plane parallel to the stream (ellipse) , as 
described above, to that in the crossflow plane (circle) . The resulting 
change in effective Reynolds number is roughly proportional to the actual 
change in streamline length (to the separation line) with angle of attack. A 
table of K a values are given below; 
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In an attempt to evaluate the independence principle for similar 

x 

boundary layer and separation conditions, the Bursnall and Loftin data is 
presented as a function of Rgff in figure 2.1. For the purpose of this 
paper, the angle of yaw has been converted to angle of attack and the drag 
coefficient defined as normal force coefficient (based on f reestream dynamic 
pressure and body diameter) in order to be more readily associated with the 
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discussion of three-dimensional bodies at angle of attack to be presented in 
subsequent sections. In order to establish the degree to which the normal 
force coefficient based on normal velocity is independent of angle of attack# 
the parameter c n /sin 2 a is plotted as a function of the effective Reynolds 
number, R e ff* for the various angles of attack, cu The results indicate that 
the critical effective Reynolds number does not vary greatly with angle of 
attack, implying that this method is reasonably accurate in predicting the 
effect of angle of attack on boundary layer transition. 

The rapid decrease of normal force with increasing Reynolds number 
combined with some degree of inaccuracy in the method of estimating R e ff 
makes it difficult to analyze, with any degree of certainty, the normal force 
data of figure 2.1 in the critical range. However, in the supercritical 
range, where the effect of Reynolds number is greatly reduced, it is clearly 
evident that the independence principle does not hold. In this range, the 

O 

results indicate that while the values of c n /sin a are about the same at the 

high angles, 75° and 90°, as would be expected by the independence principle, 
the values increase considerably as the angle of attack decreases. The 
corresponding pressure distribution for a=90°, 60°, and 30° at a value of 
R e of approximately 4.5x10^ are presented in figure 2.2. From these pressure 
distributions, it appears that as the angle of attack, a , is reduced, the 
laminar separation bubble disappears and stable lee-side vortices are formed 
as might be expected from the impulsive flow analogy (reference 2.17) for 
these finite length cylinders. These vortices produce a large normal force, 
analogous to the vortex lift of slender wings, and restrict the aft movement 
of the separation line with angle of attack. 

Subsequent to the Bursnall and Loftin study, Trimble carried out a 
similar study of oblique circular cylinders in which he extended the 
investigation to lower angles of attack. Trimble’s results, published by 
Glenny (reference 2.10), are presented in figure 2.3 in the same format as 
figure 2.1. Although, as mentioned above, data in the critical Reynolds 
number range are difficult to analyze, it is evident that the critical 
Reynolds number is somewhat lower than that of figure 1.2 and the Bursnall ^ 
and Loftin data, possibly due to a higher wind tunnel turbulence factor. 

In the supercritical range, of primary interest in this paper, the 
Trimble data is very informative. First, it will be observed that for the 
a =30° case, representing the upper Trimble case and the lower Bursnall and 
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Lof tin case* the agreement between the two experiments is exceptionally good 

(compare figures 2,1 and 2,3), At the lower angles of attack* 25° and 20°* 

Trimble 1 s data shows essentially the same values of c n /sin 2 ot* implying that 

the flow patterns are similar and the normal force value is dictated by the 

normal component of dynamic pressure, Trimble 8 s measurements made at a =5° 

are not presented in figure 2,3 since when there is a very high axial 

component of velocity* the Esch factor does not apply. However* it is of 

interest that for a=5°, those data points which are obviously in the 
• • 0 

supercritical region* the values of c n /sin a are in excellent agreement with 

those for a=20° and a=25° (see figure 7 of reference 2,10), 

A summary of the types of flow encountered in the various angle of 
attack regimes is presented in figure 2.4 for both subcritical and 
supercritical conditions. For the subcritical condition* data measured by 
Relf and Powell (reference 2,1) are used inasmuch as it appears to be the 
only data which is definitely in the subcritical range over the complete 
angle of attack range. The data of Relf and Powell has sometimes been used 
to demonstrate the applicability of the "independence principle" over the 
complete range of angle of attack. However* as pointed out by Smith et. al. 
(reference 2.13)* a careful examination of this subcritical data also 
indicates an increase in c n /sin 2 ct as the angle of attack is decreased as seen 
in figure 2.4. This variation with angle of attack appears consistent with 
the impulsive flow analogy which* from Sarpkaya 1 s data (reference 2.17)* 
predicts that for this cylinder fineness ratio of 80* the peak value of 
c n /sin a would be about 1.29 and would occur at an angle of attack of about 

6° in subcritical flow. These low angles of attack correspond to early times 
in the impulsively started flow about a cylinder normal to the flow when only 
the two starting vortices are present prior to the wake flow generated by the 
cylinder. 

For supercritical Reynolds numbers* the type of flow has a more 

2 

pronounced effect on the variation of c n /sin a with angle of attack* as shown 

by the bottom curve of figure 2.4 for a Reynolds number of 4,5x10®, At 
angles of attack between 70° and 90°, the force and pressure data of BursnaKL 
and Lof tin indicate* as might be expected* that essentially the same flow 
characteristics exist* with regard to the laminar bubble with its turbulent 
reattachment and separation leading to the characteristic unsteady and 
unstructured wake flew with no dominant vortex shedding frequency. Under 
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these conditions, the pressure forces are proportional to the normal 
component of the dynamic pressure, and c n /sin 2 a is, therefore, independent of 
angle of attack in this range. 

An analogous situation exists in the low angle of attack range where the 
combined data of Trimble and of Bursnall and Loftin indicate that in the 
angle of attack range from about 5° to about 40°, c n /sin 2 a is again 
independent of angle of attack but that its value is considerably greater 
than for the high angle of attack range. As discussed previously, this lower 
angle of attack range is characterized by a steady, well- organized vortex 
flow which restricts the rearward movement of the separation line and results 
in a large, vortex- induced increase in normal force. 

In the intermediate angle of attack range of about 40° to 70° where the 
flow might be described as "quasi-steady", there is a transition from the 
steady, well-ordered vortex flow to the unsteady, unstructured wake type 
flow. In this transition region, the parameter c n /sin a is, of course, no 
longer constant. 

Unfortunately, there is no data in the higher supercritical and 
hypercritical Reynolds number range which is the range of primary interest 
with regard to full scale aircraft and missile applications, A possible 
consequence of hypercritical Reynolds number conditions can be speculated by 
returning to Trimble® s results in figure 2.3. It is seen that for the 20° 
and 25° angles of attack, the values of c n /sin 2 a are essentially the same and 
relatively constant with Reynolds number to the highest value investigated, a 
trend which might be expected to continue due to the predominance and 
stability of the vortex flow. On the other hand, the extreme angle of attack 
case represented by the results for the cylinder normal to the stream 
(a=90°) , reproduced from figure 1.2, shows an increase in the normal force 
parameter in the middle supercritical range. In this range, the turbulent 
separation line has moved forward to a more stable position, and the wake is 
more structured and is characterized by the return of a dominant vortex 
shedding frequency (see section 3) . The drag has increased as a result of 
the forward movement of the turbulent separation and approaches the value of 
the normal force parameter at the low angles of attack. In view of the above 
effects, it is conceivable that in the hypercritical Reynolds number range, 
the normal force parameter c n /sin^t may be relatively independent of angle of 
attack over the entire range. 
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The above possibilities should be explored with the Reynolds number 
capability of the new cryogenic high Reynolds number tunnels, at least for 
three dimensional bodies, as will be discussed in section 4. 

Additional information relative to oblique circular cylinders can be 
found in references cited in section 3. 
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3, 


VORTEX SHEDDING FROM TWO-DIMENSIONAL CYLINDERS 


The wakes behind bluff bodies normal to the wind, such as the two- 
dimensional cylinders discussed in section 1, are characterized by shed 
vorticity of various frequencies which give rise to oscillating forces, 
particularly in a plane perpendicular to the wind. In aerospace engineering, 
these oscillating forces are of importance in such areas as aeroacoustic 
radiation, buffeting associated with protuberances, and oscillating loads on 
launch vehicles induced by surface winds. 

Although this paper does not deal with oscillating forces, a knowledge 
of the two-dimensional vortex-shedding frequency is helpful in understanding 
some of the steady asymmetric vortex flow characteristics of slender three- 
dimensional bodies at large angles of attack (to be discussed in section 5) . 

A brief review of the experimental studies regarding the influence of 
Reynolds number on vortex shedding beind two-dimensional bluff bodies will, 
therefore, be given in this section. 

The reference list presented at the end of the section is divided into 
subject areas and in chronological order for the convenience of those readers 
wishing to make a more thorough study. Therefore, the citations in the text 
do not follow in numerical order. 

Although the Aeolian tones produced by the action of the wind across 
wires, tree branches, etc. have been observed since antiquity, the first 
scientific research on this subject in modern times was published by Strouhal 
in 1878 (reference 3.1). Strouhal demonstrated that the Aeolian tones were 
independent of the length and tension of the wires, and depended on the 
diameter, D, and relative wind speed, V, only. He established the relation 
between these parameters which led to the now well-known “‘Strouhal number”, 

S, the non-dimensional frequency parameter given by; 


where n is the frequency of vortex shedding from one side of the cylinder. 

At essentially the same time, Raleigh concluded that the production of 
the tones was connected with the instability of the vortex sheets. He 
performed experiments which verified his suspicion that the wind induced 
vibrations of the Aeolian harp were not in the wind direction but 
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perpendicular to it and reported the results in 1879 (reference 3*2)® 

Shortly thereafter? he concluded that the Strouhal number is dependent on the 
Reynolds number (reference 3.3)? a fact that was to be clearly demonstrated 
by the research of Relf? published in 1921? and of Relf and Simmons published 
in 1924 (references 3.24 and 3.25). While their research adequately covered 
the subcritical Reynolds number range and provided a brief excursion into the 
critical range where they noted that dominant periodicity ceases? it was some 
twenty years later before significant research began in the higher Reynolds 
number regimes - research that is still ongoing but only slowly defining the 
detailed content and effects of the vortex wakes. 

3.1 Circular Cylinders in Normal Flow 

3.1.1 Vortex shedding data 

The effects of Reynolds number on the vortex shedding frequency 
associated with two-dimensional circular cylinders normal to the stream at 
low speeds are summarized in figures 3.1 and 3.2 in terms of the Strouhal 
number? versus Reynolds number based on diameter. In the interest of 
clarity and analysis? figure 3.1 presents only the results obtained from time 
histories of either lift fluctuations or surface pressure fluctuations? while 
figure 3.2 presents only those obtained from wake velocity fluctuations. The 
main features of the vortex shedding characteristics can be conveniently 
described with the aid of figure 3.1, From power spectral density analysis? 
the vortex wake character can be divided into various regimes depending upon 
the Reynolds number. In the well-known subcritical Reynolds number regime? a 
well-defined dominant periodicity exists at a Strouhal number of about 0.19? 
a value that has been obtained by many investigators from both local and wake 
measurements (the band of a large number of pressure fluctuations 
measurements made by Drescher (reference 3.31) being representative). In the 
critical Reynolds number range? where the laminar bubble with turbulent 
reattachment and separation is formed? the wake becomes narrow? and three- 
dimensional effects begin to appear in the separated shear layer due at least 
in part to the critical nature of the bubble stability. In this range? the 
dominant periodicity ceases to exist? and the fluctuations can be described 
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as wide-band random. By the time the supercritical region is reached 
(R=4xl0 5 ) , the frequencies have been observed to vary between two somewhat 
dominant values? a low value about equal to the subcritical Strouhal 
frequency, and a high value of the order of twice the subcritical frequency. 
'Phis region, defined in the figure by the data of Loiseau and Szechenyi 
(references 3,40 and 3.41) extends through the supercritical Reynolds number 
regime with only a slight decrease in the high frequency value and a slight 
increase in the low frequency value with increasing Reynolds number. 

From an analysis of local pressure spectra, Loiseau and Szechenyi 
conclude that the high frequency lift fluctuation mentioned above does not 
represent an increase in vortex shedding frequency but is the result of 
double-frequency pressures on the far leeward surface of the cylinder induced 
by the alternate shedding. They propose that as a result of the randomness 
of shedding in this Reynolds number range, a lift force at the double- 
frequency can be generated. However, in light of high frequency fluctuations 
observed in the near wake by other researchers and the implications of the 
narrowing of the wake deduced from the "universal wake Strouhal frequency" 
concept to be discussed later, it would appear to this author that high 
frequency shedding can occur in the critical and supercritical Reynolds 
number range. 

In view of the wide-band nature of the wake fluctuations in this 
Reynolds number range, the Strouhal number becomes somewhat academic for the 
two-dimensional cylinder in normal flow. However, the establishment of the 
general trend in high frequency shedding with Reynolds number may be of value 
in analyzing the cyclic shedding along three-dimensional bodies of revolution 
in the angle of attack range where there is a stable and well-organized 
wake. This three-dimensional flow is related to an impulsively started flow 
about a cylinder at early times, prior to the generation of the unstructured 
wake. It is, therefore, felt that the high level of Strouhal frequency 
observed may, in fact, be representative of the vortex shedding along three- 
dimensional bodies of revolution in the critical and supercritical Reynolds 
number ranges. This will be discussed further in section 5 in relation to ^ 
out-of-plane forces generated by three-dimensional bodies. 

Returning to the two-dimensional cylinders, Jones et. al. (reference 
3.38) found that in the upper supercritical range, from a Reynolds number of 
about 3x10® to 6x10®, where the wake has again widened, the oscillation 
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spectra is of the narrow-band random type with a single peak (at a Strouhal 
number of about 0,22.) occurring at a Reynolds number of about 3x10 s » They 
found a fairly rapid increase in Strouhal frequency as the Reynolds number 
was increased further? with a Strouhal frequency of about 0,28 being reached 
at a Reynolds number of about 6x10 s . Above this Reynolds number, they found 
a ’’quasi-- per iodic" vortex shedding condition with a value of the Strouhal 
frequency parameter of about 0.30 being reached at a Reynolds number of about 
8x10 . (The flagged symbols for the Jones et. al. data represent those 
measurement made in freon.) In this paper, the end of the supercritical 

ft 

Reynolds number range is taken to be 6x10 corresponding to the beginning of 
the "quasi-per iodic" vortex shedding region and the region above the 
supercritical is termed "hypercritical." 

Although higher Reynolds number data obtained at higher Mach numbers 
indiciates that a "plateau" may have been reached, additional data is needed 
to establish the characteristics in the hypercritical range. This indication 
of the reappearance of periodic vortex shedding in the hypercritical range is 
consistent with the earlier study by Roshoko (reference 3.15) using wake 
measurements. With regard to the reappearance of periodic vortex shedding, 
another important contribution of the study by Jones et. al. is the 
convincing evidence from both surface flow pictures and pressure measurements 
that at the highest Reynolds number, three-dimensionality in the separating 
shear layer has in general been eliminated, and the flow is again essentially 
two-dim ensi onal . 

The general characteristics observed in figure 3.1, which is based on 
force characteristics,, are essentially duplicated in figure 3.2 which is based 
on wake velocity fluctuation measurements. Each set of data in figure 3.2 
was obtained by means of a hot wire probe located at the distances, X, aft of 
the rear of the cylinder. The same regions of fluctuation characteristics 
observed in figure 3.1 are found in the wake frequency data of figure 3.2. 
Differences observed in the onset of the frequency rise in the critical range 
may be associated with tunnel turbulence or differences in model surface 
condition. With regard to surface condition and the high Reynolds number 
data obtained by Roshko (reference 3.15), the reader is reminded that the 
drag analyses (section 1.1) indicated that the sandblasted surface had 
apparently contributed a measurable drag increment because of the thin 
boundary layer in the hypercritical range. It also should be noted that the 
Belany and Sorensen supercritical data presented in figure 3.2 have been 
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limited to that obtained at or below M=0 .20 , because some of the reductions 
in the Strouhal frequency parameter found in their study at the higher 
Reynolds numbers could be the result of a compressibility effect observed in 
their drag data (see section 1.1.1). 

With regard to the high wake frequencies measured in the critical and 
supercritical range, it is recognized that controversy exists concerning the 
interpretation as in the case of the lift fluctuations described earlier. 

For example, Ericsson and Redding (reference 3.23), in questioning the high 
values of Strouhal number measured by Del any and Sorensen, speculate that 
they were associated with fluctuation in the wake formation region ahead of 
the shed vortex street because the measurements were made at positions less 
than two diameters downstream of the trailing edge. They cite the wake 
measurements of Modi and Dikshit (reference 3.22) and conclude from the 
measured vortex street geometries that the frequency could be more than twice 
as high at one diameter downstream than at 3.5 diameters downstream. It 
appears, however, that Ericcson and Redding did not account for the increase 
in vortex velocity with distance downstream found in reference 3.22, a factor 
which tends to result in a relatively constant frequency. In addition, it 
must be remembered that the Modi and Dikshit study was carried out at 
subcritical speeds where essentially all investigations have identified the 
same value of the Strouhal frequency regardless of the measuring station. 
Clarification of the vortex-street stability criterion could be of help in 
understanding the phenomena in the critical and supercritical range and this 
will be touched on in the following brief review of the universal wake 
Strouhal frequency studies. 

3.1.2 Universal wake Strouhal frequency 

While primarily of academic interest for two-dimensional circular 
cylinders normal to the flow because of the "broad-band" nature of the 
fluctuations in the wake, identification of the vortex shedding frequency in 
the critical and supercritical Reynolds number ranges may, as mentioned 
earlier, be of help in predicting the vortex shedding characteristics of 
three-dimensional bodies of revolution in the moderate and high angle of 
attack range. In these Reynolds number ranges, some additional insight might 
be derived from application of the "universal wake Strouhal frequency" 
concept, 
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Von Karxnan (reference 3.4) investigated the stability of the alternating 
vortex filaments shed from the sides of two-dimensional bluff bodies and 
found that for straight parallel paths only one arrangement was stable? that 
being the case where the distance between the two rows is 0,28 times the 
distance between successive vortices in the same row. Because of this 8 Von 
Karman street" characteristic several researchers have attempted to derive a 
universal Strouhal number, based on various wake characteristics, that would 
relate the vortex shedding frequencies of various bluff bodies. The earliest 
study appears to be that of Fage and Johansen (reference 3.7) who, by using 
measured wake widths and free stream velocity, found essentially the same 
wake Strouhal number for cylinders having rectangular (flat plate) , circular, 
and wedge shaped cross sections. The most widely used universal wake 
Strouhal number is that developed by Roshko (reference 3.11) who, recognizing 
that the Fage and Johansen formulation held only for bodies having nearly the 
same wake velocities, developed a universal wake Strouhal number based on the 


S* « Mi » 0.164, 
v s 

Von Karman vortex street stability criteria, where n is the vortex shedding 
frequency, d' is the distance between the free streamlines (shear layers) in 
the wake downstream of the coupling region and V g is the velocity along the 
free streamline at the separation point. Solving for the frequency one gets 

n _ 0.164V„ 

n — I g f 

d s 

which relates the vortex shedding frequency to the velocity and width of the 
wake. For the Strouhal frequency based on cylinder diameter, d, and for free 
stream velocity, V, 


S = 0.164 £_ Za 
d s V 


By theoretically relating d/d s and V s /V to measured drag and base pressure, 
Roshko was able to predict measured shedding frequencies reasonably well in 
the subcritical Reynolds number range. 
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From this model of the wake vortex street* one might expect that the 
large reduction of d 5 in the critical and supercritical Reynolds number range 
due to the aft movement of the separation line might result in a higher 
Strouhal frequency* S. This* of course* assumes that the reduced velocities 
at the separation point do not completely compensate for the wake width 
reduction. 

The most recent formulation of a universal Strouhal number appears to be 
that developed by Bearman (reference 3.20) which utilizes Kronauer 8 s proposed 
criterion for vortex street stability (reference 3.17). Kronauer’s criterion 
states that for a given vortex velocity the vortex street adjusts itself into 
the configuration giving minimum drag. Using this criterion* Bearman defined 
a new universal-wake Strouhal number and developed a method of predicting the 
standard body Strouhal number by use of on measured base pressure and drag. 

Inasmuch as experimental wake studies exist that appear to cast some 
doubt on the universal nature of the Karman stability criterion (see 
reference 3.20 for review)* an examination of predictions of the Strouhal 
frequency by a universal wake Strouhal frequency concept will be limited to 
several applications of the Bearman method just described. Furthermore* in 
an attempt to reduce errors due to wind tunnel interference* etc. * 
predictions were limited to those sets of data where both drag and base 
pressure were measured in the same facility. 

The results are shown in figure 3.3* and while only three sets of 
reliable base pressure and drag data for M<0.2 and covering a Reynolds number 
range were found* the prediction method does appear to substantiate that high 
levels of wake vortex shedding frequency can be expected in the Reynolds 
number range where the wake is relatively narrow. The Bearman points* 
calculated from the three tabulated values of base pressure and pressure drag 
presented in Table I of reference 3.20, appear to substantiate the trend in 
the critical Reynolds number range* while the corresponding calculations from 
the Bur snail and Loftin data (reference 2.7) extend the trend to the early 
portion of the supercritical range. In the upper supercritical range* the 
general level of the calculations based on Roshko's pressure and drag data is 
in qualitative agreement with the previously observed trend of reduced 
Strouhal frequency in this range where the wake is again widening. The 
detailed shape of the variation however is not reproduced. 

It will be noted that the base pressure and pressure drag data from the 
study by Jones et. al. (reference 3.38) have not been used here to estimate 
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the Strouhal shedding frequency. The reason for this omission is that the 
pressure measurements were made with the gap between the fixed and 
instrumented sections of the cylinder in a configuration which was found to 
seriously disturb the flow and produce unreasonably high pressure drags? 
particularly in the critical and supercritical regions. 

From the results of the estimates of the Strouhal vortex shedding 
frequency just described (figure 3.3)? it appears that the narrowing of the 
wake in the critical and supercritical regions has a greater effect than the 
reduced wake velocity? and that increases in the vortex shedding frequency- 
might be expected in this Reynolds number range. In addition? the estimated 
values are at least in qualitative agreement with the wake measurements 
presented in figure 3.2 and might reasonably be considered as additional 
evidence of the existence of high vortex shedding frequencies in the critical 
and supercritical Reynolds number ranges. However, an in-depth study of the 
applicability of the Bearman formulation to the wake characteristics in these 
high Reynolds number ranges may be required before coincidental agreement can 
be completely ruled out. 

Somewhat related to the above discussion is the study of Shaw, presented 
in a rarely cited paper (reference 3.14) in which he develops an acoustic 
theory to explain vortex shedding characteristics. Briefly, it is based on 
the assumption that the flow about the cylinder is in acoustic vibration at a 
frequency determined by the perimeter of the cylinder. He shows reasonably 
good agreement with the Bel any and Sorensen measurements for various cross- 
sectional shapes including the high Strouhal shedding frequencies in the 
supercritical Reynolds number range. However, in view of the fact that Shaw 
had to assume a harmonic mode for each case, and the fact that no independent 
critique of the basic concept appears to have been published, his results 
should not, at present, be considered as strong evidence in support of the 
high frequencies. 

3.1.3 Effect of compressibility 

Although data defining the effect of compressibility on the vortex 
shedding frequency over the Reynolds number range of interest in this paper 
is rather meager, it appears that some trends can be established as shown in 
figure 3.4. For comparison, the variation with Reynolds number for the 
incompresible case (M<0„2) is also shown. In the subcritical to 



supercritical Reynolds number range, the incompressible wake measurements of 
Bearman (reference 3,37) are shown, and in the upper supercritical range, 
those of Del any and Sorensen, In the hypercritical range, the incompressible 
lift oscillation data of Jones et. al. were selected to provide consistency 
with their compressible flow data. 

In the subcritical Reynolds number range, the compressible flow data of 
Naumann et, al. (reference 3.36) and of Murthy and Rose (reference 3.42) 
indicate the possibility of a slight reduction in Strouhal frequency with 
increasing Mach number. Once the critical Mach number is exceeded, the wake 
width tends to be fixed by shock-induced separation, and the Strouhal 
frequency remains at essentially the same value well into the Reynolds number 
range corresponding to that defined in incompressible flow as supercritical. 
This is shown by the shaded conjectural fairing in figure 3,4 based on the 
above mentioned data for a Mach number of approximately 0.45. This shock- 
induced separation precludes the wake narrowing that occurs with turbulent 
separation at subcritical Mach numbers, and therefore, the shedding frequency 
associated with shock-induced separation is essentially the same as with 
incompressible laminar separation. This, of course, results in a large 
effect of compressibility in the Reynolds number range, corresponding to the 
incompressible critical and lower supercritical ranges. In their study, 
Nauman et. al. observed that the shock waves alternated from one side of the 
cylinder to the other at the frequency of the vortex shedding. 

In the upper supercritical range, a possible trend towards reduced 
vortex shedding frequency with increasing Mach number is observed by 
comparing the Del any and Sorensen data for MC0.2 with the conjectural band 
for Mf 0.29 developed from their data and the data of James et. al. In this 
shock-free subcritical Mach number range, the apparent reduction in frequency 
with increasing Mach number is probably associated with the increase in the 
adverse pressure gradient, due to compressibility, and the resulting increase 

in the vortex wake width, as discussed in section 1.1.3. It appears that in 

\ 

the wide- band random wake region increasing Mach number reduces the upper 
frequency level. 

In the hypercritical Reynolds number range, the data point for M=0»4Q, 
obtained from Table IV of Jones et, al. , indicates that the reduction in 
Strouhal frequency with increasing Mach number in the subsonic range extends 
into this Reynolds number range. However, as in the subcritical Reynolds 
number range, the effect of Mach number is considerably less than in the 
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critical and supercritical range. It should be noted that the data point 
presented for £$=0,4 is that obtained by Jones et. al, in air with their data 
obtained in freon not being presented due to possible I! real gas i! effects for 
the high local Mach number conditions involved as discussed in section 
1 » 1 « 3 , 

Data into the supercritical Mach number range is primarily limited to 
low Reynolds numbers with the data of Murthy and Rose (see reference 3.42) 
and Knowler and Pruden (see reference 3.28) indicating that between M=0.40 
and M=0 .80 ? the Strouhal frequency remains at about a constant value. In 
this regard? it is of interest to note that the extraction of the Strouhal 
frequency number from the vortex shedding patterns of three dimensional 
bodies of revolution at high angles of attack performed by Thomson and 
Morrison (see figure 15 of reference 3.39) using the impulse flow analogy 
also indicated little variation with crossflow Mach number below M=0„8 in 
the subcritical Reynolds number range. 

In summary? it appears that for Mach numbers below about 0.80? the only 
major effects of compressibility occur in the critical and supercritical 
Reynolds number range where the narrow wake is progressively increased in 
width by the effect of compressibility on the adverse pressure gradient and 
separation point and finally by the shock-induced separation which occurs 
near the laminar separation point. For Reynolds numbers below and above this 
region? compressibility effects are minimized due to the similar location of 
the incompressible separation points and the shock-induced separation 
points. 

For Mach numbers above 0.80 very little definitive information appears 
in the literature with regard to the vortex wake characteristics for 
cylinders in normal flow. It does appear? however? that above 0.90 no 
distinct vortex shedding frequencies have teen observed in the wake (see 
reference 3.42 for example). Thomann (reference 3.34) suggests that 
supersonic flow stabilizes the separated shear layers. It might also be said 
that the fact that the wake boundaries converge rather rapidly behind 
cylinders at supersonic speeds may indicate that the lack of a distinct 
frequency might be analagous to the incompressible supercritical Reynolds 
number case with its narrow wake and no dominant vortex frequency. 

With regard to the analysis of flow about three-dimensional bodies? the 
analysis of Thomson and Morrison (see figure 15 of reference 3.39) ? described 
earlier in the subsonic discussion is of considerable interest for transonic 



and supersonic speeds* By using the impulse flow analogy to extract the 
Strouhal f requencies, they found a rapid increase of frequency with 
increasing Mach number for transonic and supersonic crossflows. Although no 
explanation was offered, it would appear that an increase in the Strouhal 
frequency number might be expected due to the reduced velocity behind the 
detached bow shock. Nevertheless, this is an area requiring further research 
which will be discussed in section 5 in connection with out-of-plane forces 
on three-dimensional bodies at high angles of attack. 

3.2 Non-Circular Cylinders in Normal Flow 

Vortex shedding frequency data for non-circular cylinders in flow normal 
to the cylinder axis is essentially limited to that of Del any and Sorensen 
(reference 3.30) and of Modi and Dikshit (reference 3.22). The former study 
covered a variety of cross-sectional shapes and in most cases included data 
through the critical range and into the low supercritical range. The latter 
study was limited to elliptical cross sections at subcritical conditions, and 
no further reference will be made to that study due to the emphasis in this 
paper on high Reynolds number. 

A summary of the effect of cross-sectional shape on the vortex shedding 
frequency parameter for a series of shapes selected from the Del any and 
Sorensen study is presented in figure 3.5. In this figure, both the Strouhal 
number and the Reynolds number, R w , are based on the maximum width, w, normal 
to the free stream velocity vector. Since Del any and Sorensen did not 
present results for the circular cylinder in the subcritical and critical 
range, the results of Bearman (reference 3.37) are utilized. The results 
indicate, as would be expected, that the effect of cross-sectional shape on 
the critical Reynolds number is consistent with that observed from the drag 
measurements (see section 1.2.1 and reference 3.30). Similarly, the effect 
of cross-sectional shape on the magnitude of the Strouhal number in the 
supercritical range follows the drag trends, with an increase in S 
accompanying a decrease in drag coefficient as might be expected from 
application of the "universal wake Strouhal frequency" concept described in 
the previous section. Unfortunately, no data appears to exist in the 
hypercritical range for non-circular cross sections. Likewise, no data 
appears to exist on the effects of compressibility for non-circular cross 
sections. 



3*3 Cylinders in Oblique Flow 


Available data on the vortex shedding frequencies of cylinders in 
oblique flow appears to be limited to the laminar separation case,, well 
below the critical Reynolds number (see references 3,44 through 3.51). Of 
the above studies? Chiu and Lienhard (reference 3.47) have provided the most 
complete coverage of the flow incidence angle range and have shewn that for 
the laminar separation case? the Strouhal vortex shedding frequency follows 
the independence principle rule which? in terms of angle of attack? a? (of 
interest for application to three-dimensional bodies of revolution) ? states 
that 


"*90 ^ si n ct ? 


where S is based on free stream velocity and body diameter and sin accounts 
for the effect of angle of attack on the free-stream velocity component 
acting in the plane normal to the body axis. 

Although no data appears to be available for the turbulent separation 
case of interest in this paper? it would be expected that the mutual 
influence of angle of attack and Reynolds number on the separation and wake 
characteristics? as in the case of the normal force, would play an important 
role in the effect of angle of attack on vortex shedding frequency. Research 
throughout the high Reynolds number range and for compressible flow is needed 
to identify such effects. This need appears to be particularly important 
with regard to vortex shedding and out-of-plane forces on three-dimensional 
bodies that will be discussed in section 5. 
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4 . FORCES ON THREE-DIMENSIONAL BODIES 


4.1 Bodies of Revolution 

4.1.1 Flow regimes 

With modern tactical aircraft and missiles encountering high angles of 
attack throughout much of their flight envelopes, the aerodynamic 
characteristics of the fuselage often play a dominant role in the flight 
behavior of the vehicle. Throughout a large portion of the angle of attack 
range, the flow about fuselage-type bodies is dominated by separation- induced 
vortex flows as illustrated by figure 4.1 (adapted from reference 4.42). 
Figure 4.1 illustrates the general trend of the boundaries separating the 
various flow regimes as a function of angle of attack and length along the 
body for ogive-cylinder types of bodies. In the low angle of attack range, 
and as long as the flow remains attached, a vortex-free flow exists over the 
length of the body, and the normal force (which is confined to the ogive) 
increases linearly with angle of attack. At an angle of attack of about 6° 
to 8°, the crossflow about the body begins to separate over the aft portion 
creating a pair of symmetrical, counter-rotating vortices fed by vortex 
sheets originating along the boundary layer separation lines on the 
cylinder. As the angle of attack increases further, the separation lines 
move toward the windward side of the body until the vortices completely 
dominate the lee-side flow. In this symmetric vortex flow regime, the 
influence of the vortices causes the normal force to increase rapidly in a 
non-linear manner similar to the vortex lift generated by highly swept 
wings. 

Above about 25°, the vortices over the aft portion become asymmetric and 
when the angle of attack typically reaches about 40°, the entire vortex flow 
is asymmetric and characterized by alternating pairs of asymmetrical 
vortices. In this steady, or at least quasi- steady, asymmetric vortex-flow 
regime, large "out- of -plane” side forces are generated, and the normal force 
no longer increases as rapidly. At even higher angles of attack, there is a 
transition into the unsteady wake-type flow characteristic of cylinders 
normal to the flow. The "out-of-plane” side forces will be discussed in 
section 5. 
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Since the flows described above are generated by boundary layer 
separation? they are influenced considerably by Reynolds number? 
compressibility? and cross-sectional shape* While a rather large mass of 
data exists for subcritical Reynolds numbers? little is available in the 
supercritical regime and essentially none in the hypercritical range. The 
purpose of this section? therefore? is to present a reference list of those 
papers which deal with the higher Reynolds numbers and to use selected 
examples to illustrate some of the influences of supercritical Reynolds 
numbers as a possible aid in guiding research programs to be carried out in 
the new high Reyonlds number wind tunnels. First? however? some of the early 
investigations of the vortex flow on bodies of revolution will be reviewed 
briefly. 

Much of the early research on bodies of revolution at angle of attack 
was related to the aerodynamic characteristics of airship hulls? and the 
"airship theory" developed by Munk (reference 4.1) ? published in 1924? is 
generally recognized as the key development in understanding the flow about 
bodies of revolution and the fundamental basis for the various slender body 
and slender wing theories subsequently developed for both attached and 
separated flows. Using a potential flow impulse analogy? he demonstrated 
that at angle of attack? the closed bodies of interest for airship hulls 
developed positive lift ahead of the maximum diameter and negative lift aft 
so that although there was no overall lift? there was an unstable pitching 
moment. However? even for the low angles of attack of interest at the time? 
some lift was developed which varied non-linearly with angle of attack? and 
Munk pointed out that this effect was associated with the fact that free 
vortices are created near the rear of the body in real viscous flow so that 
instead of a negative lift? a positive lift is produced over the aft 
portion. 

The first detailed study of the flow field associated with these 
vortices appears to be that carried out by Harrington and reported in 1935 
(reference 4.2), However? little progress was made in developing methods 
that predicted the nonlinear lift? particularly for high angles until the 
late 1940 ! s. With the growing interest in maneuverable missiles and the 
lengthening of fuselage forebodies as a result? in part? of engines being 
located in the rear of jet fighter aircraft? the high angle of attack 
characteristics of bodies were rapidly becoming more important. 
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The first successful approach to the above problem was developed by 
Allen and published in 1949 (reference 4,3) . In his approach? the 
development of the viscous crossflow along the length of the body was 
related to the time-dependent viscous flow past a two-dimensional cylinder 
started impulsively into motion with a velocity equal to Vsin . In this 
"’impulse flow analogy" distance from the nose of the body is related to the 
time from the start of the cylinder motion. 

Referring again to figure 4.1? the four flow regimes for the three- 
dimensional bodies? starting at zero angle of attack? correspond in general 
terms to the following stages of the impulsively- star ted cylinder flows (1) 
the completely attached flow? (2) the formation of two counter rotating 
symmetrical vortices? (3) the asymmetric displacement of the two vortices? 
and (4) the formation of an alternating "street 1 * of shed vortices. 

To simplify the application, the variation of the cylinder drag during 
the early stages of the motion was ignored and steady state crossflow drag 
coefficients were applied along the body. The approach? known as the Allen 
"crossflow theory? " was extensively evaluated by Allen and Perkins 
(reference 4.4)? who also studies the vortex flow fields. Some pressure 
distribution and flow field studies have been selected and are listed as 
references 4.5 through 4.12. 

The work of Allen and Perkins formed the basis for many of the semi- 
empirical methods developed since that time. While it is not the purpose of 
this paper to evaluate the various prediction methods? a list of some 
selected papers is included as references 4.13 through 4.27 for those 
interested in details of the various refinements to the basic method. 

Progress is? of course? being made with methods that take advantage of 
the rapid development of high speed computers to model the actual three- 
dimensional viscous flow. However? the current inadequacy of transitional 
and turbulent boundary- layer flow modeling and high computer costs has 
generally limited the application of these approaches to Reynolds numbers 
below those in the critical range, and no references to these studies will be 
made. 


4.1.2 Effects of angle of attack 

While a great mass of experimental data on the aerodynamic 
characteristics of bodies of revolution for a wide range of Mach numbers 
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appears in the literature* there appears to be few studies covering both a 
large angle of attack range and the subcritical* critical* and supercritical 
Reynolds number ranges. One of the most complete studies covering these 
conditions appears to be that of Hartman (reference 4.40)* and some of his 
results are presented in figure 4.2 for a high-fineness ratio* ogive-cylinder 
body at a Mach number of 0.5. The ogive has a fineness ratio of 3 and the 
cylinder a fineness ratio of 16. The results are presented in the form of 
total body normal force coefficient C N (based on free-stream conditions and 
base area) as a function of Reynolds number based on free-stream velocity and 
the diameter of the cylindrical part of the body. Due to the possibly large 
sting support interference mentioned by Hartmann* his data presented here is 
limited to an angle of attack of 60°. 

The data presented in figure 4.2 illustrates the large reduction in 
critical Reynolds number (when based on diameter) with decreasing angle of 
attack. This is associated with the increasing streamline length as 
discussed previously in connection with two-dimensional oblique cylinders in 
incompressible flow (see section 2) . However, for the case shewn* the free- 
stream Mach number is constant (M=0.5) , and since the crossflow Mach number 
varies from approximately 0.17 to 0.43 over the angle of attack range 
presented* some compressibility effects would be expected in the upper range* 
as discussed in section 1.1.3. Unfortunately* as indicated in section 1.1.3, 
two-dimensional compressibility effects in the critical Reynolds number range 
are not well enough defined to accurately apply them to the angle of attack 
case. Nevertheless* in an attempt to show at least qualitatively how changes 
in streamline length caused by changes in angle of attack can effect the 
boundary-layer separation* boundaries based on the "effective" Reynolds 
number concept discussed in section 2 are presented in figure 4.2. The 
boundary defining the beginning of the critical range has been approximated 
by R = 4x10 VlC a and that defining the beginning of the supercritical range by 

6.5xlO^/K a . Since* as mentioned above, no attempt has been made to account 
for the moderate changes in the compressibility effects due to angle of 
attack, the constants in the above approximations for the boundaries have 
been selected to provide a reasonable match with the data for the 
intermediate angle of attack of 40°. The variation of (the ratio of the 
effective Reynolds number to the diameter Reynolds number) with angle of 
attack is presented in the table in section 2.2. In general* these 
boundaries appear to be reasonably consistent with the data* although the 
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data is not sufficient to define the supercritical regime at the high angles 
of attack. In addition to the basic angle of attack effect , the expected 
effect of compressibility in delaying the beginning of the critical Reynolds 
number range appears to be in evidence at 60° angle of attack where the 
crossflow Mach number has reached 0,43, 

Another point of interest illustrated by figure 4,2 is that for a 
constant Reynolds number,, R Dr of about 3.5x10^ all three flow regimes are 
encountered - the supercritical*- critical, and subcritical - as the angle of 
attack is increased from 20° to 60°, The corresponding changes in the 
separation lines associated with the changes in effective Reynolds number 
increase the nonlinear behavior of the normal force above that associated 
with the increase due to the increase in crossflow dynamic pressure. 

To provide some indication of the variation in the upper supercritical 
Reynolds number range, data for an angle of attack of 30°, (from Foley's 
study, reference 4,34) is presented by the "flagged” symbols in figure 4.2. 1 
Since Foley's data was obtained with a fineness ratio 7 cylindrical 
afterbody, the data has been adjusted to approximate a fineness ratio of 16 
afterbody by increasing that part of the normal force acting on the afterbody 
(using slender-body theory to subtract the nose force) by the ratio of the 
fineness ratios. Despite the simplicity of the above adjustment, the 
resulting data appears to provide a reasonable extrapolation of Hartmann's 
data. For reference, the approximate boundary between the supercritical and 

fk 

lypercritical Reynolds number regimes, based on 6x10 VK^, is shewn. 

The overall results shown in figure 4.2 appear to reflect the general 
effects of Reynolds number anticipated from the two-dimensional oblique 
cylinder results described in section 2.2, although the level of the normal 
force in the upper supercritical region is somewhat higher than would be 
expected. This may be due to the fact that the method used to adjust Foley's 
data to approximate a fineness ratio of 16 may not properly adjust the 
contribution associated with the vortex flow. From the literature search, it 
appears that an experimental study covering the complete supercritical range 
and extending well into the hypercritical would be highly desirable. 


^Foley's uncorrected data is used since for this high angle condition, where 
stable vortex flow exists, it would appear that a correction based on 
asymmetries observed near zero angle of attack would not be applicable, 
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With regard to compressibility effects, it must be kept in mind that a 
more accurate estimate of the boundaries., using the Esch factor, requires 
well-defined compressible flow data for the two-dimensional cylinder case 
that would allow the effect of the normal component of Mach number to be 
included* 

Further illustrations of the significance of the flow regime encountered 
by slender bodies of revolution over the angle of attack and Reynolds number 
ranges are described with the aid of figures 4.3 and 4.4. In figure 4.3, the 
normal force data of figure 4.2 has been converted to the normal force 


coefficient on the cylindrical section based on planform area and crossflow 

8 2 

dynamic pressure, C N /sin a, and presented as a function of the effective 

H cyl 

Reynolds number, ®effr which is equal to RqK^ This provides a convenient 
way to relate the normal force to the drag of a two-dimensional cylinder in 
flow normal to its axis which is representative of the 90° angle of attack 
case. The normal force on the cylindrical section was extracted by assuming 
that the force on the expanding nose section was essentially equal to the 
mviscid value, sin2acos ^ (based on maximum cross-sectional area) and 
subtracting it from the total. The results presented for angles of attack of 
20° and 30° are seen to agree reasonably well in the critical Reynolds number 
range, indicating that the effective Reynolds number factor, K^, provides a 
reasonable means of assuring similar boundary layer conditions. The data 
from Tinling and Allen (reference 4.29) , denoted by the flagged symbols, and 
extend the a=20° data into the upper supercritical region. 

In the supercritical region, the data for the two angles of attack do 
not collapse as well, due possibly to compressibility effects. In this 
regard, the data obtained by Hartmann at the higher angles of attack are not 
presented because the normal component of Mach number for those cases were in 
the region where compressibility effects are known to be important. 

The drag coefficient variation for a circular cylinder normal to the 
flow, taken from figure 1.2, is also shown in figure 4.3. This provides an 
indication of the normal force parameter values that would be expected for the 
cylindrical portion of the three-dimensional body at a=90° where an 
unsteady wake-type flow exists. Comparison of this curve with those for the 
lower angles of attack provides an indication of the effects of the stable 
vortex flow which characterizes the moderate angle of attack range. For 
subcritical Reynolds numbers, the comparison appears to indicate that, 
consistent with the results of section 2.2 for oblique two-dimensional 
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cylinders? the effect of angle of attack is relatively small. This is also 
consistent with the impulse flow analogy as can be shown by application of 
Sarpkaya's measured? time-dependent? impulsive flow drag for a two-dimensional 
cylinder in normal flow at subcritical Reynolds numbers (reference 4,17), 

Unfortunately? no impulsive flow drag data is available for the other 
Reynolds number regimes. However? the large increase in C N ! /sin 2 a for a=20° 

a=30°? relative to a=9Q°? seen in figure 4.3 in the supercritical range is 
consistent with the oblique cylinder results described in section 2.2. Seme 
insight into the phenomena can be obtained from figure 4.4 where pressure 
distribution data? in the form of Cp/sin 2 a? is compared for a two-dimensional 
case and a three-dimensional case for supercritical Reynolds numbers. Both 
sets of data are for an effective Reynolds number of about 0.88x10 s with the 
three-dimensional body data being taken from Tinling and Allen’s study 
(reference 4.29) and the two-dimensional cylinder data taken from Bursnall 
and Loftin (reference 2.7). Excellent agreement is seen on the windward half 
of the body. However? on the leeward half of the body? it appears that the 
vortex flow restricts the aft movement of the separation line? as well as 
inducing suction forces on the leeward portion of the body. 

4.1,3 Compressibility effects 

An extensive study of the effects of compressibility over a broad range 
of Reynolds numbers has been reported by Foley in reference 4.34. The 
effects of Mach number and Reynolds number on normal force coefficient for 
the body at an angle of attack of 30°? taken from the above study? are 
presented in figure 4.5. The free stream Mach numbers range from 0.40 to 
2.0? while the normal component of Mach number? M sin a? varies from 0.20 to 
1.0. As in figure 4.2? the normal force coefficient is based? in the 
conventional manner for bodies alone, on the base area. Thus? because of the 
lower fineness ratio and moderate angle of attack (relative to figure 4.2) ? 
the coefficient level is greatly reduced and a different scale is selected to 
better illustrate the compressibility effects. 

The general trends of the compressibility effects on this three- 
dimensional body at 30° angle of attack are the same as those observed in the 
drag data for two-dimensional cylinders normal to the flow (see section 
1.1.3). However? the magnitude of the forces illustrate the effect of the 
stable vortex flow that occurs on the three-dimensional bodies. Sizable 
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compressibility effects are observed? even in the subcri tical crossflow Mach 
number range. Also? the large increase in normal force that occurs above a 
crossflow Mach number of 0,4? accompanied by a decrease in Reynolds number 
effect? is consistent with the shock induced separation conditions observed 
in the two-dimensional flow. 

Again it appears that the research needs to be extended into the 
hypercritical Reynolds number range, A useful contribution would be to 
develop accurate Euler solutions for the inviscid shock induced separation 
case that would provide an asymptote for the very high Reynolds number 
conditions. Such solutions would provide useful information regarding the 
need for additional experiments to extend the shock boundary layer 
interaction results to higher Reynolds numbers. 

Additional papers containing the effects of Reynolds number and Mach 
number on the aerodynamic characteristics of bodies of revolution can be 
found in the reference lists of both this section and section 5. 

4.2 Bodies of Non-Circular Cross Section 

4.2.1 Normal force characteristics 

The cross sections of fuselages often depart from circles? and an 
example of the important effects of these shapes on the normal force is shown 
in figure 4.6, which is based on data obtained on fuselages similar to the 
space shuttle orbiter and reported by Brownson et. al. in reference 4.31. 

The data is for an angle of attack of 60° and covers a range of windward 
corner radii for cross sections of modified rectangular shapes. The Mach 
number has been limited to values less than 0.2 to minimize compressibility 
effects. The normal force coefficient is based on the projected plan 
area of the fuselages? and the Reynolds number is based on the maximum width? 
W? which is more closely related to the boundary layer development ahead of 
the corner than is the body depth. The radius parameter? r/ w ? covers the 
range from 0 to .086, In the subcri tical range? it will be noted that only 
small increases in radius produce large normal force reductions relative to 
the sharp corner case. The corresponding decrease in critical Reynolds 
number? analagous to that encountered in two-dimensional flow (see section 
1.2.1)? emphasizes the importance of high Reynolds- number testing when 
investigating the implications of windward corner shape. The effects of 
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corner shape are of particular importance with regard to the local normal 
force and its impact on pitching moment. 

Some effects of angle of attack and Mach number for the same shuttle 
type fuselage are presented in figures 4,7 and 4.8, which show data taken 
from the summary paper by Jorgensen and Brownson (reference 4,35), For these 
comparisons, the fuselage having r/^,063 has teen selected. Figure 4.7 

presents the data for a Mach number of 0,3 in the form of C^/sin^a as a 
function of an effective Reynolds number. In selecting the effective 
Reynolds number factor, it appeared probable that due to the reduced lateral 
perturbations, the streamline lengths would increase more rapidly with 
decreasing angle of attack for this flat bottom fuselage than for the 
circular fuselages. Therefore, instead of using K a (section 2.2) , the factor 
1/sina was used and as shown in figure 4.7, was found to collapse the data 
for angles of attack from 35° to 50° as well as could be expected using this 
simplified approach to the three-dimensional boundary layer development. For 
this angle of attack range, the normal component of Mach number, M sina is 
less than 0.23, so that the data is assumed to be reasonably free of 
compressibility effects. For the angle of attack of 70° where the value of 
M sina has reached 0.28, there appears to be a consistent effect that may be 
attributed to compressibility promoting a more forward separation line in the 
vicinity of the small-radius corner with its steep adverse pressure 
gradient. 

Regarding the effective Reynolds number factor selection, it should be 
noted that Esch (reference 4.38) selected l/sin 2 a and found a reasonably good 
correlation over the complete angle of attack range. However, the present 
author believes that the larger factor required in Esch 8 s analysis is due to 
the probable mixture of subcritical and supercritical Mach numbers 
encountered in the high angle of attack range for the data he used. 

To further illustrate the effects of compressibility, figure 4.8 
presents the normal force parameter, C N i /sin 2 a as a function of I^/sina for 

free- st ream Mach numbers of 0,2, 0.6, and 0,9, for angles of attack of 40° 
and 70°. For convenience in relating compressibility effects for the two 
angles of attack, the normal component of Mach number, M sina, is also 
shown. The most prominent effect shown is the much greater influence of 
compressibility at the angle of attack of 40°, where vortex flow would be 
expected, than at 70°, where turbulent wake flow is expected. This 
difference in compressibility effect is particularly noteworthy in view of 
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the considerably lower levels of the normal component of Mach number for the 
40° angle of attack case. The fact that the vortex-induced force begins to 
increase at a relatively low value of M sina is presumed to be associated 
with the high velocity at the separation line for the small corner radius. 

For both angles of attack, the reduction in Reynolds number effects is 
observed for Mach numbers where shock induced separation would be expected. 

In general, the primary importance of the fuselage normal force 
variation with Reynolds number is, of course, the contribution it makes to 
the aircraft pitching moment. An example of the effect of Reynolds number on 
the pitching moment of a forward-swept wing aircraft configuration at low 
speeds and an angle of attack of 80°, taken from unpublished data obtained at 
NASA 1 s Ames Research Center is presented in figure 4,9. In view of the long 
fuselage forebody, and the fact that the remainder of the fuselage is 
"blanketed” by the canard and wing, it is to be expected that the major 
effect of Reynolds number would be associated with the fuselage forebody. It 
is seen that as the Reynolds number increases, there is a large negative 
shift of the pitching moment consistent with the expected reduction in 
forebody normal force as the crossflow attaches further around the 
forebody. Such effects are important with regard to stall-spin 
characteristics, and some related data on this configuration can be found in 
reference 4.61. 

The results presented in figure 4.9 are for cross sections with rather 
generous radii on the bottom corners? for smaller radii the critical Reynolds 
number would, of course, be higher, as indicated in figure 4.6. 

The above results tend to emphasize the importance of simulating both 
the Mach number and the Reynolds number of the flight conditions in question, 
even for relatively low Mach number conditions, for these types of body cross 
sections. The new cryogenic wind tunnels offer important increased 
capability to accomplish this well into the hypercritical Reynolds number 
range where no data exists at present. 

4.2.2 Side force characteristics 

One of the most pronounced effects of non-circular bodies is the side 
forces they develop, particularly when the sideslip is combined with the high 
angles of attack encountered in maneuvers. For tactical aircraft, the major 
effects are usually confined to the fuselage forebody ahead of the wing and 
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the resulting side forces often produce yawing moments of magnitudes such that 
the complete aircraft can become statically unstable. 

Section 1.2.2 demonstrated the importance of Reynolds number on the side 
forces developed on various two-dimensional non-circular cylinders with axes 
normal to the flow and at roll angles representing various combinations of 
angle of attack and sideslip. The implications of these side forces with 
regard to three-dimensional bodies were briefly discussed. Also discussed 
were some of the early aircraft configuration studies of non-circular 
fuselages which demonstrated the need for two- and three-dimensional 
research, and the reader is referred to section 1.2.2 for this background and 
associated references. 

This section briefly describes some Reynolds number studies carried out 
on three-dimensional bodies in the high angle of attack range. Extensive 
studies on three-dimensional bodies having non-circular cross sections have 
been performed by Clarkson et. al. (reference 4.43) over an angle of attack 
range from 55° to 90°, and a portion of their results are presented in figure 
4.10. The side force coefficient, Cy‘ , acting on the forward portion of the 
fuselage is presented as a function of for a constant flow incidence angle 

in the body crossflow plane of <f>=+10° (obtained by static roll angle) for 
various angles of attack. The angle 4> represents the resultant flow angle in 
the crossflow plane for various combinations of angle of attack and 
sideslip. 

The fuselage cross section was a square with rounded corners having non- 
dimensional radii of ^=0.245 which is identical to one of the two- 
dimensional cylinders discussed in section 1.2,2, For angles of attack of 75 
and 90° the resulting side force coefficients, based on side area, are quite 
similar to the two-dimensional results presented in figure 1.11, and 
illustrate the change in direction of the side force between subcritical and 
supercritical Reynolds number conditions as the flow attaches around the left 
front corner (see section A-A of figure 4.10), Inasmuch as a positive value 
of corresponds to a positive sideslip angle, it is seen that the nose 
portion would contribute a restoring moment at subcritical Reynolds numbers 
and diverging moment at supercritical Reynolds numbers. For configurations 
where the wing obstructs the crossflow aft of the nose, such effects as 
shown in figure 4,10 can have a profound effect on directional stability, and 
even more important, at these high angles of attack, on damping of the flat 
spin mode (to be shown later) . 
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For the angles of attack of 45° and 60° , large Reynolds number effects 
still occur but are seen to have somewhat opposite trends. It would appear 
that this change is associated with the change from a disorganized wake type 
flow to the organized vortex type flow discussed in connection with circular 
bodies. For example, at the lower angles of attack, a strong vortex would be 
expected to form and be fed by the separation at the left front corner, in 
the cross-sectional view of figure 4*10, which produces suction forces 
accompanied by circulation changes such that the resulting force is opposite 
that at the higher angles of attack. At higher Reynolds numbers where the 
flow remains attached further around the corner, the vortex strength is 
reduced and a reduction in side force occurs. 

An example of the Reynolds number sensitivity of space shuttle orbiter 
type bodies at high angles of attack combined with sideslip is presented in 
figure 4.11, from the study by Brownson et. al. (reference 4.31). Again, the 
side force coefficient, Cy* , is presented as a function of the Reynolds 

number based on the maximum width of the body. The angle of attack is 60°, 
the flow incidence angle, <f> , is 10°, and the non-dimensional corner radii 
r/w, is 0.086. Also shown is the corresponding data from the two-dimensional 
cylinder study, reviewed in section 1.2.2, for the square cylinder having non- 
dimensional corner radii r/w of 0.080. The two-dimensional side force has 
been adjusted to the condition of ot=60° by reducing its magnitude by sin^a to 
account for the reduced crossflow dynamic pressure. Since the effective 
Reynolds number for the a=60° case is higher than for the a =90° case, the 
Reynolds numbers for the two-dimensional data have been reduced for purposes 
of comparison. It should be noted, however, that instead of the Ka term 
presented in section 2.2 for circular bodies, the factor sina has been used 
as discussed in section 4.2.1. Although such corrections are very 
approximate, and changes in the lee-side flow can occur that are not 
accounted for, the two sets of data demonstrate the same trend of side force 
with Reynolds number. 

The major consequence of the above side force characteristics is 
usually, of course, its effect on the vehicle's directional stability. An 
example of such an effect of Reynolds number on the directional stability, 

C n , is presented in figure 4.12 for the complete configuration of the space 

shuttle orbiter. This data is for an angle of attack of approximately 20° 
and was taken from the paper by Bornemann et. al. (reference 4,41). The data 
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points represent static directional stability measurements made in various 
wind tunnels with several models which resulted in a variety of stability 
levels. However*- when the data was plotted as a function of the Reynolds 
numbers associated with the various size models and tunnel pressures* the 
trend illustrated in figure 4,12 emerged. The results indicate that while 
the configuration is stable at low Reynolds numbers, it becomes unstable at 
high Reynolds numbers. This again illustrates the importance, particularly 
for non-circular forebodies, of testing at sufficiently high Reynolds numbers 
to assure separation characteristics similar to the full scale vehicle, or of 
developing accurate transition-control techniques. It will be noted that the 
change in directional stability is consistent with the side force results 
presented in the previous figure, for a body alone having a cross section 
similar to the shuttle oribter nose. Although the body alone data was for a 
much higher angle of attack, it appears possible that with the upwash from 
the wing increasing the local angle of attack at the nose section, side force 
variations similar to those of the previous figure could be developed on the 
nose, and that they would provide the observed trend in directional 
stability. The long moment arm from the center of gravity amplifies, of 
course, the effect of the side force developed on the nose. It should be 
pointed out that these results include any forces induced on the aft portion 
of the vehicle by the vortex flows generated by the forebody separations. 

The new high Reynolds number facilities offer an opportunity to carry 
out a systematic study of this and other viscous dominated flows on the space 
shuttle orbiter and other configurations at, or near, full scale conditions. 

4.2.3 Autorotation characteristics 

Another consequence of the long fuselage forebodies ahead of the wing, 
which are characteristic of modern jet fighter-type aircraft, is the major 
role they can play in regards to the aircraft* s spin characteristics. Prior 
to the evolution of these long forebodies, the fuelage’s influence in 
spinning was usually small. However, in the mid-1950 9 s, several aerodynamic 
research programs began to demonstrate not only the strong influence of the 
longer forebodies but the critical effects of their cross-sectional shape and 
of Reynolds number. The early portions of this research were influenced 
considerably by the discovery in the Langley spin tunnel that the model of 
the F8U aircraft experienced a flat fast spin. These programs included* in 
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addition to the spin tunnel studies, one-degree- of- freedom autorotation tests 
of the F8U in the flat spin mode over a range of Reynolds numbers and a 
fundamental study of the crossflow characteristics of a wide range of two- 
dimensional non-circular cylinders over a wide Reynolds number range, These 
latter two programs demonstrated the strong influence of cross-section shape 
and Reynolds number and the results relative to the F8U are briefly 
summarized with the aid of figure 4,13, From their two-dimensional cylinder 
research, Polhamus et. al, (reference 4,49) established a probable 
autorotation boundary for rectangular and square cross sections with rounded 
corners as a function of the Reynolds number (based on the corner radius and 
the ratio of the corner radius to the width of the cross section) as shown in 
figure 4.13. In the same time period, Clarkson (reference 4.48) was carrying 
out the F8U autorotation tests in which he established a Reynolds number that 
separated autorotation from damping charcteristics for the fuselage forebody 
and his results, converted to the parameters described above, are also shown 
in figure 4.13. The lower Reynolds number of the boundary for the F8U cross 
section relative to the modified rectangular cross sections is probably due 
to a less adverse pressure gradient expected for the somewhat elliptical 
shape as well as three-dimensional relief. Both studies predicted the 
importance of Reynolds number which was substantiated by comparisons of the 
F8U spin tunnel results by Kliner et. al. (reference 4.44) and the full scale 
aircraft flight results by Alley (reference 4.45). These results also shown 
in figure 4.13, indicated that, in fact, the airplane was not susceptible to 
the flat spin observed in the spin tunnel. 

As a result of the above studies, the current practice in the spin tunnel 
tests is to develop a forebody flow control device for each configuration to 
simulate high Reynolds number damping characteristics. This is accomplished 
through the aid of static tests carried out over a Reynolds number range in a 
conventional pressure tunnel as has been described by Chambers et, al. 
(reference 4.51). Some of the results of these static tests can be found in 
references 4.52 through 4.55, Although this approach has been fairly 
successful, uncertainties associated with the lack of rotary motion in the 
Reynolds number tests and possible incompa labilities between the flow control 
devices required for the yaw and pitch modes are causes for some concern. 

Recently, Clarkson et, al. have established a general research program 
utilizing a one-degree-of-f reedom rig similar to the earlier F8U rig, but 
with angle of attack variation capability. This program includes force and 
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moment data? detailed pressure distribution measurements? and surface flew 
visualization for a series of body shapes over a Reynolds number range. This 
technique should offer? in addition to the basic research results? an 
improved method of flow control development for use in improving simulation 
in low Reynolds number spin tunnels. An example of the side force 
(indicative of the damping or propelling moment) measured on this rig? over a 
range of Reynolds numbers is shown in figure 4.14 for a fuselage having a 
square cross section with non-dimensional corner radii? r/w ? of 0.245 (see 
reference 4.56). Also shown is the side force measured on a two-dimensional 
cylinder of the same cross section for a flow incidence angle? <j> ? of 10° 

(from figure 1.11). The corresponding helix angle for the rotational tests 
was 7.6° and the strong similarity of the results is apparent. Both types of 
tests indicate a shift from a prospin condition at low Reynolds numbers to an 
antispin condition at higher Reynolds numbers as the flow attaches around the 
left front corner in the cross-section sketch. A detailed description of the 
flow is presented in section 1.2.2. Additional results from the spin rig can 
be found in references 4.56 through 4.58 

As pointed out in section 1.2.2? cross-sectional shapes such as the 
modified squares might be expected to encounter additional Reynolds number 
effects in the higher ranges as the flow eventually closes more completely 
around the aft corners. Therefore? it would be advantageous to carry out 
tests similar to those described above in the new high Reynolds number 
tunnels and obtain actual full scale conditions. 
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5 


OUT-OF-PLANE FORCE 


5*1 Background 

5*1*1 Early studies 

As the angle of attack of a body of revolution increases, the crossflow 
separation induced vortex pairs eventually become asymmetric even for the 
symmetric body at zero sideslip. This can result in a large out-of-plane 
side force producing a sudden and sometimes uncontrollable yawing moment. 

Some of the early studies are listed as references 5,1 through 5.12. 

The first photographs of these asymmetric vortices to be published were 
those of Allen and Perkins in 1951 (reference 5.1). They pointed out that on 
long missiles at high angles of attack, "the discharge of a vortex street 
should promote asymmetry of forces on the tail surfaces and manifest itself 
as a tendency to unexpected and erratic rolling ..." The anticipated effect 
on long missiles with stabilizing fins was soon confirmed by wind tunnel 
tests as reported by Mead in early 1952 (reference 5.2) and in 1953 by Gowen 
(reference 5.3) , who reported on vertical tail buffeting caused by asymmetric 
vortex flow. 

With regard to forces and moments developed on the body itself, Allen 
and Perkins also suggested that "the aperiodic changes in the discharged 
street of vortices might induce, for a pitched body, undesirable forces and 
moments in yaw. 81 In the latter part of 1952, Cooper et» al, published the 
results of a study at M=l„59 which demonstated the out-of-plane forces and 
moments anticipated by Allen and Perkins and, in addition, indicated the 
cyclic variation in local side force along the body length (reference 5.4). 

In 1953, Letko (reference 5,5) published the first low speed measurements of 
out-of-plane forces and moments on slender bodies at high angles of attack 
and demonstrated the effectiveness of a ring or other roughness at the nose 
in greatly reducing the undesirable yawing moments. Also in 1953, Perkins 
and Dennis (reference 5.6) established the angle of attack at which the 
asymmetric and unsteady vortex flow began to decrease with increasing nose 
fineness ratio, and then established a boundary relating these two 
parameters, A study of the yawing moment at zero sideslip was also made in 
connection with spin stabilized rockets by Dunn in 1954 (reference 5.9). 
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Some early work in England was reported by Fink in 1953 (reference 5,11) 
and by Mai tby and Peckham in 1956 (reference 5,12), During this same period? 
studies established that? although the asymmetric vortex patterns were still 
present? the large out-of-plane forces and moments encountered at low speeds 
were essentially eliminated at supersonic speeds. Representative of this 
research are the papers published by Mien and Perkins (reference 5,1) ? Gowen 
and Perkins (reference 5.8)? and Gapcynski (reference 5.10), 

Thus? the general features of the asymmetric vortex flow? the resulting 
out-of-plane forces? the cyclic variation of local side force, and the 
importance of nose shape and Mach number had been uncovered by the mid 
1950 8 s. However detailed research, the development of prediction methods, 
and the search for practical methods of reducing the magnitude of the forces 
lagged during the next fifteen years because the aircraft and missiles of 
that time were not? in general, operating in the high angle of attack range. 
However, in the late 1960*s and early 1970 s s, interest was renewed as new 
fighter aircraft designs began to penetrate the high angle of attack range 
during maneuvers. At the high angles of attack involved the large yawing 
moment which can result from the side force at zero sideslip can be a major 
contributor to the so-called "nose slice" or "yaw off" - a divergence in yaw 
resulting in an incipient departure. Providing adequate control power for 
these high angle of attack asymmetric flow conditions is a major problem for 
modern tactical aircraft that can encounter angles of attack approaching 60° 
in combat maneuvers. This and other stability and control problems 
encountered in the design of modern fighter aircraft have recently been 
reviewed by Skow and Erickson (reference 5.20). Examples of review papers 
dealing with the fluid dynamic aspects of the asymmetric vortex flows and 
resulting out-of-plane forces are listed as references 5.13 through 5.22 with 
the recent paper by Hunt (reference 5.21) providing a summary of the overall 
subject. 

with regard to basic data papers, the reference list has been restricted 
primarily to those that include studies of Reynolds number and/or Mach number 
effects in the subsonic and transonic range (consistent with the subject of 
the present paper) . The large mass of low speed, subcritical Reynolds number 
data, as well as supersonic data, can be located within the reference lists 
of the cited papers. 

Because of the many recent papers dealing with various aspects of this 
flow, and because of the overall summary by Hunt, the discussion in the 
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present paper will be limited primarily to certain effects of Reynolds number 
and a discussion of those areas where higher Reynolds number data might be 
needed. 

Although the resulting yawing moments are generally more important to 
vehicle controllability , the side forces and their distributions have been 
selected as a more convenient means of discussing the general phenomena 
associated with asymmetric vortex flow. 

5.1.2 Some general characteristics 

A typical variation of the out-of-plane force with angle of attack for 
an ogive-cylinder configuration is shown in figure 5,1. At a moderate angle 
of attack,, vortex flow develops but is initially symmetrical. At some higher 
angle of attack,, the vortices become asymmetrical with the asymmetry either 
being initiated at the rear of a relatively long afterbody in accordance with 
the impulse flow analogy or at the apex of a slender, pointed nose. The side 
force resulting from the vortex asymmetry increases with angle of attack due 
primarily to the increase in the normal component of dynamic pressure. 
However, the corresponding increase in normal velocity increases the vortex 
shedding rate, resulting in a second pair of asymmetric vortices of opposite 
asymmetry, and the overall side force is reduced. Eventually, as the angle 
of attack increases further, the steady vortex flow transitions into the 
unsteady wake-type flow associated with cylinders normal to the flow, and the 
time-averaged side force returns to zero. 

For aircraft-type bodies, it appears that since the wing generally 
blankets the flow over much of the fuselage, the onset of vortex asymmetry is 
associated with the apex of the nose, as mentioned above for short bodies. 
Keener and Chapman (reference 5.34) have proposed that since asymmetric 
vortex flows and the resulting rolling moments have been observed in tests of 
slender sharp-edge wings where the separation lines are fixed, asymmetry in 
separation lines on slender pointed bodies of revolution are not necessarily 
the cause of vortex asymmetry ; instead it can result from a "hydrodynamic 
instability" of the vortices themselves. Woolard (reference 5.38) arrived at 
a similar conclusion by showing that the onset angles for slender wings and 
cones can be analytically related by a slender body mapping procedure. More 
recently. Dyer et, ai, (reference 5.39) have shown, with an inviscid model of 
the flew about a slender cone, that asymmetrical vortices can form with the 
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separation lines disposed symmetrically. These and other papers related to 
hydrodynamic instability are listed as in references 5,29 through 5,40, 

Figure 5,2 illustrates some of the evidence of "hydrodynamic 
instability" as a cause for vortex asymmetry. On the left are results of 
tests of slender delta wings that shew the boundary between symmetric and 
asymmetric vortex flow in terms of a/6 * the ratio of angle of attack to semi- 
apex angles versus the semi-apex angle. On the right are the results of the 
inviscid theory for cones in terms of a/6 versus the location of the assumed 
symmetrical separation lines. Experimental data for cones indicate the 
boundary occurs at a value of | of about 2.0* which is in reasonable 
agreement with the calculations of Dyer et, al. and about one- half of the 
value for the delta wings. 

From these results and available experimental data* it appears that the 
onset angle of attack may not be critically dependent upon Reynolds number. 
However* as will be shown, the characteristics beyond onset are highly 
Reynolds number dependent. 

5.2 Ogive- Cylinder at Constant a 

The recent study by Lamont (reference 5,57) produced a rather complete 

set of data with regard to the detailed effects of Reynolds number on the out- 

of-plane force for an ogive-cylinder covering the subcritical* critical* and 

supercritical regimes. While only a portion of his study has been published 

to date* that which has been published provides considerable insight into the 

Reynolds number efects of interest in this paper. The variation of the 

maximum out-of-plane force, C Y * with Reynolds number* R D * for an angle 

x o max 

of attack of 40° from Lamont's study is presented in figure 5.3. The model 

was a body of revolution having an ogive nose of 3.5 diameters in length and 

a cylindrical afterbody of 4.0 diameters in length of which 6 diameters from 

the apex were pressure instrumented. The model was rolled during the tests 

to assure that the maximum side force was obtained. 

It will be noted that in preparing figure 5.3* this author utilized the 

actual sign of the side force rather than using the absolute value* IC Y I 

I x o max I 

presented by Lamont. The parameter* |C V I* is often useful when 

I 1 o maxi 

comparing results from different models or different runs since the direction 
of the vortex asymmetry and the resulting side force sign are generally 
established by random disturbances in the wind tunnel or slight asymmetries 
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in nose shape. This method of presentation? however? obscures the changes in 
sign of the overall side force that are caused by changes in Reynolds 
number. These changes are sometimes encountered in the critical range and 
are believed to be associated with the change in vortex shedding frequency. 

The actual signs used in figure 5.3 were deduced from the longitudinal 
distributions of local side force. 

In figure 5.3, it will be noted that C Y decreases with increasing 

o max 

Reynolds number, reaching a value of zero at about R D =7xlO^ and then changes 
sign. With further increases in Reynolds number, the sign again reverses, 
and at the highest supercritical Reynolds number tested, %f=4xl0 6 , the out-of- 
plane force returns to a value similar to the subcritical value. While much 
of the reduction in the critical range can be attributed to the well-known 
reduction in vortex strength as the separation line moves aft, the change of 
sign appears to be associated with the change in period of the cyclic 
variation of local side force along the length of the body. This is 
illustrated in figure 5.3 by the small inset plots showing the measured 
longitudinal distributions of the local side force for three Reynolds 
numbers. It will be noted that near the peak negative value of out-of-plane 
force, the period of the cyclic variation has been reduced to about one-half 
of its subcritical value and the negative loop establishes the overall sign 
of the side force. In the upper supercritical Reynolds number range, the 
period returns to a value similar to the subcritical case, and a return to a 
positive value of side force coefficient results. The cyclic variation of 
side force along with body is, of course, associated with the alternating 
asymmetric vortex pairs resulting from periodic vortex shedding, and it can 
be shown by the impulsive flow analogy that the vortex spacing varies 
inversely as S tana/ where S is the Strouhal frequency for a two-dimensional 
cylinder normal to the flow (see reference 5.61) having the same separation 
line location. A1 though the analogy has been successfully demonstrated in 
the subcritical range (see reference 5.50, for example), no such analysis 
covering a wide range of Reynolds numbers in which the various flow regimes 
are encountered has been found in the literature. However, the detailed data 
set of Lamont's offers an opportunity to explore this aspect and the result 
of a preliminary attempt is illustrated in figure 5.4. 

In figure 5.4, Lamont's local side force distributions along the body 
are presented for several Reynolds numbers on the left portion of the 
figure. The Reynolds numbers were selected such that they could be 
74 



reasonably characterized by the laminar? transitional? or fully turbulent 
separation boundaries identified by Lamount (see figure 2? reference 5*57) . 
The corresponding two-dimensional Strouhal vortex shedding parameter? S? was 
then estimated from figure 3,2 of this paper. While the choice of a value 
for the shedding parameter was straightforward for the laminar and high 
Reynolds number turbulent cases? the values to use for the transitional and 
low Reynolds number turbulent cases were not so apparent because of the 
variation of the boundary-layer state along the length of the body. 
Consequently? for the purpose of this paper? the shedding frequency 
parameters for the latter two cases were chosen as those values appropriate 
for the longitudinal station at which the first peak in the side force 
distribution occurs. Consequently? using Bearman 9 s data? a subcritical value 
of S equal to 0.20 was used for the laminar separation case? a value in the 
center of the critical range equal to 0.32 for the transitional case? and a 
value in the supercritical range of 0.46 for the low Reynolds number fully 
turbulent case. A value of 0.26 based on Roshko's data was used for the high 
Reynolds number fully turbulent case. It should be noted that Lamont's 
boundary layer states were used rather than R D since it is difficult to 
relate the effective Reynolds number on the ogive to the two-dimensional 
case. 

Using the values of S shown in the table? the ratio of Cy/c n is plotted 

versus (S)tana on the right portion of figure 5.4. The use of Cy/c n is an 
attempt to reduce the effect of vortex strength variations with Reynolds 
number. While the amplitudes do not collapse? it is interesting to note that 
the periods? in terms of ^ (S)tana are essentially the same. From this it 
appears that the three-dimensional body data may support the suggestion made 
in section 3.1? that while the two-dimensional vortex shedding frequency 
content in the critical and most of the supercritical range is of the "wide- 
band random" type? the upper frequency level characterizes shedding from 
three-dimensional bodies at angles of attack where steady asymmetric vortex 
flow exists. This conclusion? however? must be considered as tentative since 
the above use of the impulse flow analogy is strictly applicable only for 
constant diameter bodies and the method of using the two-dimensional data of 
figure 3.2 could be questioned. 

To further investigate the above observations? it would appear desirable 
to make pressure measurements on a high fineness ratio ogive cylinder in 
order to define the cyclic variations of the out-of-plane side force over the 
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constant diameter section in the supercritical and hypercritical Reynolds 
number ranges* The new high Reynolds number cryogenic tunnel offers the 
opportunity to accomplish this over a range of subsonic and transonic Mach 
numbers. 

5.3 Isolated Ogive Forebodies 

5.3.1 General characteristics 

The above discussion dealt with the effect of Reynolds number on C v at 

o 

a constant angle of attack. However, additional insight with regard to the 
effects of Reynolds number can be derived by examining how the first and 

second peaks in the cyclic variation of C Y change with angle of attack. 

*0 

Inasmuch as the sense of the vortex asymmetry is random and dependent on 
disturbances associated with model imperfections, tunnel flow, etc., a sign 
convention has been adopted as an aid in clarifying trends. In this 
convention, all signs are reversed, when necessary, to provide a positive 
sign for the first peak in the subcritical Reynolds number range. An example 
of the effect of Reynolds number on these peak values is presented in figure 
5.5 for a fineness ratio 5.0 tangent ogive forebody alone at tfc=0.25. The 
data illustrate the basic variation with angle of attack for a subcritical 
and supercritical Reynolds number as obtained from the data of Keener et. 
al. (figure 6 of reference 5.52) and defines the first and second peaks. 

For the subcritical Reynolds number case shown, as angle of attack increases 
it appears that two opposing effects occur that play dominant roles in the 
variation of the out-of-plane force with angle of attack. One effect is the 
increase in crossflow dynamic pressure, which is proportional to sin 2 a, while 
the other is an increase in vortex shedding frequency which, according to the 
impulse flow analogy, is proportional to tan . The latter effect manifests 
itself, of course, as a reduced longitudinal spacing between the shed 
vortices (see reference 5.60, for example). In the angle of attack range 
where asymmetrical vortices first start to develop, the increasing crossflow 
dynamic pressure dominates and the out-of-plane force increases with angle of 
attack. However, the increase in vortex shedding frequency with increasing 
angle of attack compresses the cyclic variation of the out-of-plane force 
along the isolated forebody length and the resulting growth in the reversed 
force loop over the aft, larger diameter, portion of the forebody eventually 
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reverses the trend with angle of attack and for the case shown, the overall 
out-of-plane force is essentially zero when an angle of attack of about 65° 
is reached, and at higher angles a two-dimensional wake type of flow 
develops. 

In the critical and low supercritical Reynolds number range, the reduced 
vortex strength associated with the leeward movement of the separation lines 
tends to reduce the out-of-plane force. However, as illustrated by the 
supercritical case of figure 5,5, an additional effect combines with those 
described in the previous paragraph to further change the character of the 
variation with angle of attack. This effect is assumed to be an increase in 
vortex shedding frequency due to supercritical separation as discussed above 
in connection with figure 5.3. As a result of this additional increase in 
vortex shedding frequency and the associated further compression of the 
cyclic force variation, the counteracting reversed force loop becomes 
prominent at a lower angle of attack and above about 35° completely dominates 
with the overall out-of-plane force changing sign and the maximum force now 
occurring at the peak of the resulting second force loop. It will be noted 
that this effect of Reynolds number is consistent with Lamont's results for a 
constant angle of attack of 40° shown in figure 5.3. 

To illustrate further the effect of Reynolds number on the peak out-of- 
plane forces, figure 5.6 presents the variation of the first and second peaks 
in the force variation as defined by figure 5.5. The variation of these two 
peaks with Reynolds number has been extracted from the angle of attack data 
of reference 5.52. Also shewn, as a vertical band, is an indication of the 
critical Reynolds number as defined from the normal force data of reference 
5.52. The results illustrate the decrease in the magnitude of the first peak 
with increasing Reynolds number and illustrate the reversed sign and 
magnitude increase for the second peak as Reynolds number is increased in the 
supercritical region. The large reduction of the first peak with increasing 
Reynolds number and the opposite sign and magnitude increase for the second 
peak with increasing Reynolds number appear consistent with the above 
discussion of the interacting effects of crossflow dynamic pressure and 
vortex shedding frequency. It is recognized, of course, that the details of 
these trends are modulated by the complicated interactions between the vortex 
flew, boundary layer state, and separation line variations over expanding 
bodies. For details of the boundary layer changes, the reader is referred to 
Lamont’s recent paper (reference 5,57). 
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In view of the large variation in the magnitude of the second peak in 
the supercritical Reynolds number range and the fact that the two-dimensional 
studies described in section 3.1.1 suggest additional changes in the vortex 
shedding characteristics in the hypercritical range , research to higher 
Reynolds numbers appear desirable. 

5.3.2 Effect of fineness ratio 

The effect of fineness ratio on the C v variation with Reynolds number 

*o 

is illustrated in figure 5.7 which is based on the data of Keener et, al. 
(references 5.51 and 5.52) , The data for the fineness ratio of 5.0 is the same 
as that presented in figure 5.6 except that? in the interest of clarity, the 
data for the first peak at supercritical Reynolds numbers is eliminated since 
the peak absolute values occurred at the second peak. In the subcritical 
range reducing the fineness ratio to 3.5 appears to have only minor effects 
implying that as a result of the low frequency of vortex shedding at 
subcritical Reynolds numbers? no appreciable reversed side force loop was 
developed on either body. Other effects such as differences in body side 
area? crossflow dynamic pressure, etc. can also contribute? however? and 
further study is required before a complete understanding can be attained. 

In the supercritical Reynolds number range? there is a well-established trend 
with the reduced fineness ratio body producing considerably less side force. 
Aside from the fact that the reversed side-force loop on the low fineness 
ratio body begins its growth at a higher angle of attack, possibly limiting 
the growth, it is again not clear how much of a role the aforementioned 
factors play in this effect, 

5.3.3 Effect of cross section 

Although in general, the forebodies of aircraft fuselages are not 
axisymmetric or even composed of sheared circular cross sections? research 
relative to the out-of-plane force has been devoted almost entirely to the 
axisymmetric case. While some research has been performed? for example that 
of Krouse (reference 5.41), Keener et. al. (reference 5.52), and Skow and 
Erickson (reference 5,20), it appears that it is limited primarily to force 
data and that only Keener et. al. have covered a significant Reynolds number 
range for non-axisymmetric bodies. 
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In their study? Keener et, al, investigated the effects of modifying 
their circular tangent ogive forebodies to elliptical cross-sectional shapes 
for a case with major axis vertical (narrow side-on) and one with major axis 
horizontal (broad side-on) . Their results have been summarized by Hunt in 
Figure 72 of reference 5*21 in terms of the maximum side force coefficient as 
a function of Reynolds number. Expecting the side force? for the same maximum 
width? to be greatest for the ellipse with the "narrow side-on" Hunt 
described as an anomaly the implication that it developed? in general? the 
lowest side force. However? there are several factors? not discssed by Hunt? 
which tend to shed additional light on the experimental results. 

First? as can be seen from the two-dimensional data of Del any and 
Sorensen (reference 1.45) and Polhamus et. al. (reference 1.52)? the critical 
Reynolds number for the ellipse with the major axis vertical (particularly 
when based on maximum width as in Hunt’s figure) is considerably lower than 
that with the major axis horizontal. This shift? due primarily to the 
differences in magnitude of the adverse pressure gradients? greatly amplifies 
the differences in side force at a given Reynolds number since subcritical 
and critical boundary layer conditions are being compared in Hunt’s figure 
72. Second? the strengths of the shed vortices are greater? for the same 
body width for the ellipse with the "broad side-on" due to the higher 
velocities at the separation lines. The greater vortex strength will? of 
course? tend to offset the effect of the reduced side area. Other factors 
include the difference in planform fineness ratio which can result in changes 
in the angle of attack for the onset of vortex asymmetry and the associated 
differences in "crossflow" dynamic pressure at the angle of attack for 
maximum side force. 

In an attempt to clarify the effects of cross section? particularly with 
regard to the Reynolds number effects of interest in this paper? the data of 
Keener et, al. are analyzed in figures 5,8 and 5.9, Several factors were 
considered in making the comparisons. The tangent ogive forebody having an 
elliptical "cross section" reported by Keener et. al. was designed such that 
with the major axis horizontal? it has the same planform as the fineness 
ratio 3,5 tangent-ogive forebody of circular cross section? and with the 
major axis vertical? it has the same planform as the fineness ratio 5,0 
circular cross-section body. Therefore? in attempting to illustrate the 
effect of cross-sectional shape? the comparisons are made herein for the same 

planform fineness ratios and with the side force coefficient based on the 
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area of a circle of diameter equal to the ellipse planform width at the 
base* Also* to provide some degree of consistency in relation to the length 
of streamline run between the stagnation line and the beginning of the 
adverse pressure gradient for the different shapes* the Reynolds numbers are 
based on the maximum body dimension (in plane of major axis) at the base in 
all cases. Finally* following the approach used in the previous section* it 
is indicated whether the maximum side force shown occurred at the first or 
second peak in the variation with the angle of attack. 

In view of the difficulties involved in establishing that the maximum 
side force for each condition is indeed being encountered, the lack of 
pressure distribution measurements* questions associated with cyclic 
variations resulting from vortex shedding frequencies* etc.* the data 
presentation and discussion which follows must be considered a tentative 
attempt to establish some trends with cross section. 

Figure 5.8 presents the comparison between the circular and narrow side- 
on elliptical cross sections for the tangent ogive forebodies of fineness 
ratio 5.0. For the subcritical case where the separation lines would be 
expected to be near the maximum body width location for both forebodies, the 
narrow side-on elliptical shape generated the largest side force coefficient 
with the increase being approximately equal to the side area ratio. At 
supercritical Reynolds numbers* the maximum side force coefficients occur at 
the second peak because of the increased vortex shedding frequency* and the 
magnitude for the narrow side- on elliptical cross section is now considerably 
lower than for the circular cross section. This reduction is presumed to be 
associated with the rather large reduction in vortex strength due to the more 
complete flow closure for the more "streamlined" shape in the turbulent 
separation region* which more than offsets the larger side area. The benefit 
of a more "streamlined" cross section was suggested during the early studies 
of the problem in 1951 by Allen and Perkins (reference 5.1) as a possible 
solution to the out-of-plane force problem for configurations for which 
inclined flight occurs essentially in one plane. However* for modern* highly 
maneuverable aircraft, flight at high angles of attack is not confined to one 
plane* and the reduced directional stability associated with the elliptical 
nose with major axis vertical is* of course* highly undesirable. 

Although the narrow side-on ellipse is of little interest for aircraft 
forebodies* the results of figure 5.8 illustrate the potential peril in 

attempting to accurately establish the effects of cross section and to select 
an 



forebody shape if the data is confined to subcritical Reynolds numbers* It 
is also possible that the trends obtained in the supercritical range are not 
representative of full-scale conditions because there could possibly be 
further changes in vortex shedding frequency in the hypercritical Reynolds 
number range. These possible effects need to be investigated in the new high 
Reynolds number wind tunnel up to full-scale conditions. 

A comparison of the circular and "broad side-on" elliptical cross 
sections for the tangent ogive forebodies of fineness ratio 3.5 is presented 
in figure 5.9. To be noted is the higher critical Reynolds number for the 
elliptic cross section associated primarily with the more severe adverse 
pressure gradient, and the reduction in maximum side force at both 
subcritical and supercritical Reynolds numbers relative to the circular cross 
section. Also to be noted is the reduced magnitude of the maximum side force 
at supercritical Reynolds numbers. Again, the need for full-scale Reynolds 
number data is illustrated, and such tests should include the combined 
effects of Reynolds number and Mach number to assure adequate simulation of 
the vortex shedding frequency and the separation line locations. 

Comparing figures 5.8 and 5.9 in the supercritical range, the much 
larger reduction in maximum side force for the "narrow side-on" elliptical 
cross section than for the "broad side-on" elliptical cross section when 
compared with their corresponding circular cross section forebodies is 
striking. This benefit of the more "streamlined" shape, which has been 
discussed earlier, has sometimes been interpreted as inconsistent with the 
benefits in reduced out-of-plane force obtained with the "shark nose" 
configuration of the recent Tiger Shark version of the F-5 aircraft (see 
reference 5.20). It would appear, however, that the primary source of the 
"shark nose" benefit is derived from the very large forebody apex angle, an 
early established trend with regard to delaying the onset of vortex 
asymmetry. The major benefit of the accompanying flat elliptical cross 
section, aside from "area rule" considerations, is presumed to be an 
improvement in directional stability. 

5.4 Compressibility Effects 

As mentioned earlier, studies in the early 1950 8 s demonstrated that 
while asymmetric vortex flow also appears at supersonic speeds, the resulting 
out-of-plane forces and moments were essentially eliminated. More recently, 
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studies have been made in the high subsonic and transonic range, of 

particular interest in this paper, as well as the supersonic range. Briefly, 

these studies established that depending upon the body nose shape and 

fineness ratio, the side force tends to vanish at free stream Mach numbers as 

low as 0.8, and as expected, that the reduction was primarily related to the 

crossflow Mach number, M since. This is illustrated by figure 5.10 which has 

been prepared utilizing the subsonic and transonic data on an ogive-cylinder 

configuration reported by Pick (reference 5.42). Presented is the maximum 

out- of -plane force C v , as a function of both freest ream Mach number, M, 

*o max 

and its normal component, M sina, where the value of a is that corresponding 

to C v at each Mach number. The results indicate a rapid decrease in the 

x o max 

out-of-plane force occurring when the normal component of Mach number reaches 
approximately 0.4 which is close to the critical Mach number for two- 
dimensional circular cylinders. The force rapidly decreases and is 
essentially eliminated for normal components of Mach number, M sina, in the 
range from 0.5 to 0.6. Similar trends have been observed by other 
investigators and Hunt's survey paper (reference 5.21) presents results from 
several of these studies. Of the various phenomena which can contribute to 
this reduction in side force (see reference 5.42), it is likely that the 
reduced sensitivity of the separation lines to wake effects for the case of 
shock- induced separation is the primary cause. 

Inasmuch as the angle of attack at which the onset of vortex asymmetry 
occurs is strongly related to the forebody nose angle, it follows that the 
onset of compressibility effects in terms of freest ream Mach number should be 
dependent upon forebody nose angle for a given angle of attack. Such a 
relationship has been demonstrated by Keener et, al. (reference 5.48) from 
tests of various forebody shapes. For example, their results indicated that 
the reduction in maximum side force due to compressibility effects was 
delayed to a considerably higher freest ream Mach number as the nose fineness 
ratio was increased from 3.5 to 5.0. The fineness ratio 5.0 nose, with its 
smaller nose angle and resultant lower angle of attack for onset of vortex 
asymmetry, reaches a higher freest ream Mach number before the crossflow Mach 
number, M sina, reaches the onset of important compressibility effects. From 
this, they suggest that since the onset of vortex asymmetry has been shewn to 
occur at an angle of attack approximately equal to twice the nose semi-apex 
angle, the Mach number at which the side force has reduced to near zero can 
be approximated from the expression M sin(2S) *0.5. 
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For the free- stream Mach numbers where strong compressibility effects 
were encountered? it appears that the Reynolds numbers were subcritical in 
the above study* There also seems to be little? if any? data in the 
literature at combined supercritical Mach numbers and supercritical Reynolds 
numbers for cases in which the nose was rotated to assure the attainment of 
maximum side force* Since the angle of attack for onset of vortex asymmetry 
appears to be relatively insensitive to Reynolds number, it might follow that 
the free- stream Mach number corresponding to supercritical crossflow Mach 
number will also be relatively insensitive to Reynolds number. However, it 
should be kept in mind that this comment relative to Reynolds number effects 
is based on forebody alone results, and that with a sizable afterbody 
present, such phenomena as the cyclic variation in vortex asymmetry along the 
body length could produce interacting effects of Mach number and Reynolds 
number, particularly in the critical Reynolds number range* Again, the new 
high Reynolds number wind tunnel facilities should provide insight into these 
possible interactions. 
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6. SOME RECOMMENDED RESEARCH 


Throughout the paper some of the areas where additional research would 
be useful have been pointed out, particularly with regard to extensions to 
higher Reynolds numbers. As a possible aid in defining research projects for 
the National Transonic Facility on bodies at high angles of attack, this 
section itemizes some of the more pertinent research areas mentioned in the 
paper. To assist in relating the recommended research items to particular 
NTF research models, they are grouped under "bodies of revolution" and 
"bodies of non-circular cross section". Also, under each body type, the 
subjects addressed follow the order of the main body of the paper. 

Priorities, therefore, are not implied by the order of the research 
recommendation. It is anticipated that, in general, possible adoption of any 
of the recommendations and assignment of priorities will be influenced by 
models currently planned for the NTF that could be adapted and future models 
will, of course, be influenced by aircraft technology needs at that time. 

In addition to the three-dimensional research, there are several areas 
where fundamental research using two-dimensional models would be beneficial. 
However, in view of the initial emphasis on three-dimensional testing in the 
NTF, no definitive recommendations will be made here with regard to two- 
dimensional research except to point out that, in the absence of wall 
turntables and balances, some fundamental research could be performed with 
cylinders at fixed flow incidence angles and with surface pressure orifices. 
Some fruitful areas of two-dimensional research can be selected from the 
sections on cylinders in the main body of the paper. 

Bodies of Revolution 

Since the forces normal to the axis of revolution for any combination of 
angle of attack and sideslip can be resolved from the normal force and out-of- 
plane side force variation with angle of attack, only those two forces will 
be discussed. 

.Normal Force 

Although the major effects of Reynolds number on normal force occur in 
the critical range, the review indicated evidence of rather large effects in 
90 



the supercritical range at subcritical Mach numbers that are not well defined 
at present. Essentially* no data is available in the hypercritical Reynolds 
number range. It is* therefore* recommended that an ogive- cylinder bocjy-of~ 
revolution model having a fineness ratio of approxmately 10 be tested over a 
wide Reynolds number range with emphasis on the supercritical and 
hypercritical regimes for angles of attack up to at least 45°. The study 
should include compressibility effects for normal components of Mach number 
up to at least 0.60 at the high angles of attack. 

An aid to this experimental research would be the development of 
accurate Euler solutions for the inviscid, shock-induced separation case to 
provide an indication of the asymptote for the very high Reynolds number 
conditions prior to eventual high Reynolds number Navier-Stokes solutions. 

Qut-Qf-Plane Force 

The analysis of two-dimensional and three-dimensional vortex shedding 
evidence* as related by the impulsive flow analogy, leads to a tentative 
conclusion that in the supercritical Reynolds number range, the vortex 
shedding rate for the three-dimensional ogive-cylinder is related to the 
upper level of the wide-band random frequency spectra for the two-dimensional 
circular cylinder normal to the flow. For this case, the Strouhal number is 
more than twice the subcritical value. It is, therefore, recommended that an 
ogive-cylinder pressure model of at least 10 diameters in length be tested 
over a wide Reynolds number range in an attempt to more completely define the 
changes in the period of the cyclic variation of the out-of-plane side force 
along the body Reynolds number well into the hypercritical range. Also, 
large effects of Mach number are expected in the supercritical Reynolds 
number range and normal components of freest ream Mach number of at least 0.5 
should be investigated in the 40° to 45° angle of attack range. The 
compressibility study should be extended into the hypercritical Reynolds 
number range if model loads will allow. 

Bodies of Non -circular Cross Section 

Bodies of non-circular cross section are highly sensitive to Reynolds 
number effects, and the data and analysis presented in the main body of this 
paper indicates many areas where research at higher Reynolds numbers would be 
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valuable. Although a few suggestions for extending the research are made in 
this section? final selections for models and test programs should? of 
course? be made by an examination of the many two-dimensional and three- 
dimensional flow examples given in the body of the paper in the light of 
cross-sectional shapes of practical interest at the time, 

■ Mo -pn a L. fq jxjs 

As in the case of the bodies of revolution? little data exists in the 
upper supercritical Reynolds number range and beyond. An even greater 
shortage of data exists with regard to compressibility effects which are 
particularly important due to the higher flow acceleration encountered around 
non-circular bodies. While large effects of Reynolds number on the normal 
force produced by non-circular bodies are illustrated in the body of the 
paper? it is believed that? in general? the shape selected for Reynolds 
number and compressibility studies in the NTF should probably be based on the 
more critical effects associated with asymmetric flows discussed in the 
following sections. 

■Si de force. D ue to j aifleslip 

One of the most pronounced Reynolds number dependent flows encountered 
by bodies at high angles of attack is that associated with non-circular 
bodies in sideslip. The analysis of both two-dimensional and three- 
dimensional data illustrated that changes in the sign of the side force with 
Reynolds number can occur and data for the space shuttle orbiter demonstrated 
that such force changes on the forebody are of sufficient magnitude to cause 
a directional instability as the Reynolds number is increased. Of particular 
interest from a research standpoint is the square cross section having 
rounded corners for which the two-dimensional analysis indicated the 
possibility of additional changes in sign of the sideforce? at Reynolds 
number beyond those for which data is available? due to asymmetric closure of 
the flow around the leeward corners as influenced by vortex flow generated at 
the windward corners. It is? therefore? recommended that studies at high 
Reynolds numbers on bodies of non-circular cross section include? at least? a 
square shape with rounded corners and that consideration be given to forebody 
tests on a complete aircraft model as well as a relatively high fineness- 
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ratio body alone. Tests of the shuttle orbiter model in the NTF should, of 
course, include directional stability studies over a wide Reynolds number and 
Mach number range. 

Out-Qf-Plane Force 

Except for certain elliptical cases, there appears to be little 
information on the Reynolds number dependence on the side force at zero 
sideslip generated by bodies of non-circular cross sections at high angles of 
attack. Those non-circular cross sections which develop large Reynolds 
number dependent side forces due to sideslip would be expected to be 
particularly sensitive to the asymmetric flow produced at angles of attack 
where the lee-side vortices become asymmetric. In addition, the review of 
the two-dimensional cylinder data indicated a large effect of cross-sectional 
shape on the Strouhal frequency in the critical Reynolds number range which 
might be expected to influence the vortex shedding and, therefore, the period 
of the cyclic variation of the side force along the body length. In view of 
these anticipated Reynolds number effects and their importance with regard to 
high angle of attack flying qualities, it would appear desirable to 
investigate at least one non-circular cross section body as part of the early 
research in the NTF in order to better define a more complete research 
program. 
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SYMBOLS 


A P 

c 

c d 


:d p 


c n 

C n 

C N 


C • 


*-N 

3 

s 

% 


c 1 

y 


c x 


L o max 


D 

h 

L 


speed of sound 
base area 
planform area 

wing mean- aerodynamic chord 

drag coefficient of two-dimensional cylinders, lift for unit 
length/qW 

pressure drag coefficients 

lift coefficient of two-dimensional circular cylinders, lift for 
unit length/qD 

pitching-moment coefficient, pitching moment /qf-^c 

normal-force coefficient of two-dimensional circular cylinders, 
normal force per unit length/qD 

yawing-moment coefficient, yawing moment/qS w c 

normal-force coefficient of three-dimensional bodies, normal 
force/qA^ 

normal-force coefficient of three-dimensional bodies, normal 
force/qAp 

directional stability parameter (rate of change of yawing 
moment with sideslip) 

pressure coefficient 

base-pressure coefficient 

sideforce coefficient in body axis, side force per unit 
length/ qh 

sideforce coefficient at zero sideslip (out-of-plane force 
coefficient), side force /qA b 

maximum out-of-plane force coefficient 

cylinder diameter 

body height or spoiler projection 

body length 
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M 

n 

q 

r 

r d 

R eff 

R h 

R r 

s 

V 

w 

X 

a 

8 

6 

<i> 

0 

9 s 

n 

p 

y 

OJ 


distance to forebody geometric centroid 
Mach number, V/a 

frequency of vortex shedding from one side of body 
dynamic pressure, 

corner radius, or sand roughness height 
Reynolds number based on maximum diameter, 

Reynolds number based on an effective streamline length (see 
section 2.2.2) 

Reynolds number based on body height, 

Reynolds number based on corner radius, 

Reynolds number based on body width, — 

Strouhal number, jjp or ^ 
wing area 

freest ream velocity 

maximum width of body (D for circle) 

distance from body nose, or distance behind cylinder 

angle of attack 

angle of sideslip 

azimuthal angle defining spoiler location, also semi-apex angle 

flow incidence in plane normal to axis of cylinder 

azimuthal angle (0=0 is on windward side) 

azimuthal angle defining location of separation line 

drag proportionality factor ( In^ te^sp^^Srag" } 

free- stream air density 

free-stream air viscosity 

rolling angular velocity 
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A 


Figure 1.1 


Identification of flow regimes to be discussed relative 
to circular cylinders. 



Figure 1.2 Drag coefficient for smooth circular cylinders. 
M = 0.15. 


97 



(r/D)x.10 5 





Achenbach 
(Ref. 1.21) 

) 



J 

1 



>22 

i 


i ^ i 

jRoshho 
Ref. 1.19) 

[ 1 

Szechenyi (Ref. 1.24) 

1 1 


Smooth (fig. 1.2) 


Figure 1.3 Effect of distributed roughness on the drag 
coefficient of circular cylinders. M = 0.1 


M <0.15 

( fig. 1.2 ) 


Effect of Mach number on drag coefficient of smooth 
circular cylinders. 

































invlacid theory 


2 Train 4 > 


Figure 1.11 Influence of corner radius on the side force coefficient 
developed on modified square cylinders (reference 1.51 
and 1.52). 
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Figure 1.12 


Effect of leeward shape on the side force coefficient of 
modified rectangular cylinders. 
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Figure 1.13 Effect of windward corner radius on side force spike 
developed by modified triangular cylinders. 



Lockwood (Rof. 1.S3) 


Clr. Cyllndor 


Figure 1.14 


Influence of slight contour distortion on side force 
coefficient of circular cylinder. 
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(fig. 1.2) a = 90 
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Figure 2,3 Effect of angle of attack on normal force coefficient 
parameter for circular cylinder. 
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Figure 2,4 Variation of circular cylinder normal force coefficient 
parameter, with angle of attack for subcritical and 
supercritical conditions. 
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Figure 3,5 Influence of cross section on Strouhal number. 
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Figure 4.1 Approximate flow regimes for slender bodies of 
revolution. Adapted from reference 4.42. 



Figure 4.2 Effect of angle of attack on the normal force 
coefficient for a slender body of revolution. 

M = 0.5. Data from Hartmann (reference 4.40). 
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Figure 4.3 Indication of vortex contribution to normal force coefficient 
parameter. 
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Figure 4.4 Effect of vortex flow on pressure distribution. 
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Figure 4.5 Effect of Mach number on normal force coefficient for 
slender body of revolution. Data from Foley (reference 
4.34). 
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Figure 4.6 Effect of corner radius on normal force coefficient for 
shuttle orbi ter- type fuselage. Data from Brownson, 
et. al . (reference 4.31) . 
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Figure 4.7 Variation of normal force parameter with an effective 
Reynolds number for shuttle orbiter-type fuselage, 
r/w = .063. Data from Jorgensen and Brownson 
(reference 4.35). 
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Figure 4.8 Variation of normal force coefficient parameter with 
Mach number for shuttle orbiter-type fuselage, 
r/w = .063. Data from Jorgensen and Brownson 
(reference 4. 35) . 
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Pitching moment coefficient for forward swept wing 
configuration. Data from Whipple et. al . (reference 
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Figure 4.10 Effect of angle of attack on side force coefficient 
for body with modified square cross section. Data 
from Clarkson et. al . (reference 4.43). 
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Figure 4.11 Side force coefficient for shuttle orbiter type fuselage 
Data from Brownson et. al . (reference 4.31). 
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Figure 4.12 Directional stability variation for shuttle orbiter 
Data from Bornemann and Surber (reference 4.41). 
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Figure 4.13 Autorotation boundaries for some non-circular cross- 
section bodies. 



Figure 4,14 Forebody side force coefficient produced by spin motion. 
Data from Clarkson et.al. (reference 4.56). 
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Figure 5.1 Typical variation of out-of-plane force and vortex patterns 
with angle of attack. 
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Figure 5.2 Examples of hydrodynamic instability of vortex pairs shed 
from symmetrical separation lines. 
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Variation of local and integrated maximum out-of-plane 
force coefficient for an ogive-cylinder. Data from 
Lamont (reference 5.57). 
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Application of impulsive flow analogy to the out-of-plane 
force variation along body length. Basic data from 
Lamont (reference 5.57). 
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Illustration of first and second peaks in out-of-plane 
force variation with angle of attack. Data from Keener 
et. al . (reference 5.52). 
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Variation of first and second out-of-plane force peaks 
Data from Keener et. al . (reference 5.52). 
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Effect of fineness ratio on out-of-plane force 
characteristics. Data from Keener et. al . (references 
5.51 and 5.52). 


Effect of cross section on out-of-plane force 
characteristics for fineness ratio 5 forebody. Data 
from Keener et.al. (reference 5.52). 
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Figure 5.9 Effect of cross section on out-of-plane force 

characteristics for fineness ratio 3.5 forebody. 
Data from Keener et. al . (reference 5.51). 
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Figure 5.10 Variation of maximum out-of-plane force with Mach 
number at subcritical Reynolds number. Data from 
Pick (reference 5.42). 
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